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ABSTRACT 


Orthogonal  frequency  division  multiplexing  (OFDM)  is  being  successfully  used 
in  numerous  applications.  It  was  chosen  for  IEEE  802.11a  wireless  local  area  network 
(WLAN)  standard,  and  it  is  being  considered  for  the  fourth- generation  mobile  communi¬ 
cation  systems.  Along  with  its  many  attractive  features,  OFDM  has  some  principal  draw¬ 
backs.  Sensitivity  to  frequency  errors  is  the  most  dominant  of  these  drawbacks.  In  this 
thesis,  the  frequency  offset  and  phase  noise  effects  on  OFDM  based  communication  sys¬ 
tems  are  investigated  under  a  variety  of  channel  conditions  covering  both  indoor  and  out¬ 
door  environments.  The  simulation  performance  results  of  the  OFDM  system  for  these 
channels  are  presented. 
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EXECUTIVE  SUMMARY 


Due  to  the  fact  that  the  need  to  send  more  information  increases  each  day,  the 
military  is  also  in  need  of  solutions  to  make  more  information  transmission  to  the  battle¬ 
field  possible.  The  military  is  looking  wireless  ADHOC  networks  as  a  solution  to  this 
problem  and  orthogonal  frequency  division  multiplexing  (OFDM)  is  the  technique  used 
in  ADHOC  networks. 

OFDM  is  being  successfully  used  in  numerous  applications.  It  was  chosen  for 
IEEE  802.1  la  wireless  local  area  network  (WLAN)  standard,  and  it  is  being  considered 
for  the  fourth- generation  mobile  communication  systems.  Along  with  its  many  attractive 
features,  OFDM  has  some  principal  drawbacks.  Sensitivity  to  frequency  errors  is  the 
most  dominant  of  these  drawbacks. 

The  objective  of  this  thesis  was  to  investigate  the  frequency  offset  and  phase 
noise  effects  on  OFDM  based  communication  systems  to  determine  the  ways  by  which 
we  can  use  this  technique  in  battlefield  applications  more  efficiently.  A  comprehensive 
background  of  OFDM  based  communication  systems  is  presented,  and  the  sub -carrier 
orthogonality  requirement,  inverse  fast  Fourier  transform  and  fast  Fourier  transform 
processes  are  discussed.  Theoretical  analysis  of  frequency  errors  and  their  effects  on  the 
overall  system  performance  are  examined  based  on  the  work  reported  in  the  literature. 

Matlab  simulations  of  an  OFDM  communication  system  were  performed  for  the 
following  channels:  ideal,  additive  white  Gaussian  noise  (AWGN),  and  mobile.  The  mo¬ 
bile  channel  selected  includes  multipath  and  Doppler  effects  in  AWGN.  The  effects  of 
frequency  offset  and  phase  noise  were  studied  for  all  channel  models  by  observing  the 
performance  deviation  from  the  AWGN  case  as  a  function  of  increasing  frequency  offset 
and  phase  noise.  The  bit  error  rate  (BER)  performance  plots  of  the  OFDM  system  af¬ 
fected  by  frequency  errors  under  AWGN  and  mobile  channel  conditions  were  obtained 
through  simulations. 


xvii 


It  was  observed  that  the  OFDM  system’s  performance  degraded  as  a  power  of  the 
relative  frequency  offset.  A  relative  frequency  offset  of  0.1%  required  a  bit  energy  to 
noise  power  spectral  density  ratio  ( Eb  /Nn)  value  of  1  dB  more  than  the  AWGN  case 

without  offset,  and  a  frequency  offset  of  0.2%  required  2.8  dB  more  Eh/Na  for  the  same 
BER. 

In  the  phase  noise  case,  it  was  observed  that  the  OFDM  performance  is  degraded 
linearly  as  we  increased  the  normalised  bandwidth,  as  was  suggested  in  the  literature.  A 
normalized  bandwidth  of  0.07%,  representing  the  amount  of  phase  noise,  required  an 
Eb/N0  value  of  1  dB  more  than  the  AWGN  case  without  phase  error,  and  each  addi¬ 
tional  0.07%  increment  in  the  normalized  bandwidth  resulted  in  1  dB  more  Eh  /N0  re¬ 
quirement. 

The  effects  of  mobility  were  also  studied.  The  mobility  is  simulated  by  including 
the  Doppler  frequency  shift  and  multipath  effects.  The  simulations  suggested  that  the  sig¬ 
nal  to  noise  ratio  degradation  due  to  the  Doppler  shift,  most  of  the  time,  is  negligible.  The 
degrading  effects  of  multipath  are  observed  by  comparing  the  system  performance  under 
mobile  channels  with  that  under  AWGN  conditions. 


I.  INTRODUCTION 


Orthogonal  frequency  division  multiplexing  (OFDM)  is  being  successfully  used 
in  numerous  applications,  such  as  European  digital  audio  broadcasting  and  digital  video 
broadcasting  systems  [1,  2].  In  1999,  the  IEEE  802.11a  working  group  chose  OFDM  for 
their  5- GHz  band  wireless  local  area  network  (WLAN)  standard,  which  supports  a  vari¬ 
able  bit  rate  from  6  to  54  Mbps.  OFDM  was  also  one  of  the  promising  candidates  for  the 
European  third- generation  personal  communications  system  (universal  mobile  telecom¬ 
munication  system).  However,  it  was  not  approved  since  the  code  division  multiple  ac¬ 
cess  (CDMA)-based  proposals  received  more  support.  OFDM  is  now  being  considered 
for  the  fourth- generation  mobile  communication  systems  [3].  Therefore,  OFDM’s  per¬ 
formance  in  mobile  and  fading  environments  is  the  topic  of  many  current  studies.  These 
studies  are  also  helpful  to  determine  the  ways  by  which  we  can  use  this  technique  in  bat¬ 
tlefield  applications  more  efficiently. 

A.  OBJECTIVE 

Although  OFDM  seems  to  be  a  solution  to  keep  up  with  the  demand  of  increasing 
data  rates,  it  has  some  drawbacks.  Sensitivity  to  frequency  errors  is  the  most  significant 
of  these  drawbacks.  The  main  objective  of  this  thesis  was  to  investigate  and  document  the 
effects  of  frequency  offset  and  phase  noise  on  the  performance  of  OFDM  based  digital 
communications  under  different  channel  conditions. 

A  step-by-step  approach  was  adopted  in  order  to  achieve  the  objective  of  this  the¬ 
sis.  The  first  step  is  to  provide  a  basic  background  on  the  principles  of  OFDM.  The  rea¬ 
sons  for  the  frequency  errors  and  a  theoretical  analysis  of  these  effects  on  OFDM  systems 
are  documented.  Finally,  a  communication  system  utilizing  OFDM  was  simulated  in 
Matlab  and  the  system’s  behavior  under  different  channel  conditions  was  observed. 

To  be  able  to  observe  the  system  behavior,  the  simulation  results  for  different 
channel  models  are  presented  in  graphical  form.  Next,  the  simulation  results  obtained  in 
this  work  are  compared  to  the  simulation  results  reported  in  related  studies. 
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B.  RELATED  RESEARCH 

Due  to  its  many  attractive  features,  OFDM  has  received  much  attention  in  the 
wireless  communications  research  communities.  Numerous  studies  have  been  performed 
to  investigate  its  performance  and  applicability  to  many  different  environments.  Below 
are  some  of  the  many  studies  conducted  concerning  the  effect  of  frequency  errors  on 
OFDM  systems. 

Many  of  the  published  studies  about  the  frequency  errors  use  two  main  refer¬ 
ences.  The  first  is  the  study  of  Pollet  el  al.  on  sensitivity  of  OFDM  systems  to  frequency 
offset  and  Wiener  phase  noise  [4],  and  the  second  is  the  study  of  Moose  on  a  technique 
for  OFDM  frequency  offset  correction  [5]. 

Other  related  studies  include  the  study  of  Armada  on  the  phase  noise  and  sub- 
carrier  spacing  effects  on  OFDM  systems’  performance  [6],  the  study  of  Xiong  about  the 
effect  of  Doppler  frequency  shift,  frequency  offset,  and  phase  noise  on  OFDM  receivers’ 
performance  [7]  and  the  study  of  Zhao  el  al.  on  the  sensitivity  of  OFDM  systems  to  Dop¬ 
pler  shift  and  carrier  frequency  errors  [8]. 

C.  ORGANIZATION  OF  THE  THESIS 

Chapter  II  starts  with  an  introduction  to  OFDM  and  proceeds  to  present  the  fun¬ 
damentals  of  OFDM.  The  generation  of  an  OFDM  symbol  is  also  explained  using  a  step- 
by-step  approach.  Chapter  HI  discusses  the  frequency  errors  including  the  frequency  off¬ 
set  and  the  phase  noise.  A  simple  OFDM  signal  representation  is  adopted  to  easily  repre¬ 
sent  the  frequency  offset  and  Doppler  frequency  shift  effect  in  analytical  expressions. 

The  reasons  for  frequency  errors  are  explained. 

Chapter  IV  develops  the  Matlab  simulations  to  study  the  system  performance  un¬ 
der  different  environments.  It  also  presents  the  results.  Chapter  V  provides  a  summary  of 
the  thesis,  the  conclusions  and  suggestions  for  future  work. 

Appendix  A  outlines  the  Matlab  code  for  the  transmitter,  the  channel  and  the  re¬ 
ceiver.  The  function  of  each  m-file  and  its  relationship  to  other  functions  are  discussed. 

At  the  end  of  each  section,  the  hierarchy  of  the  m- files  is  illustrated  in  a  schematic  form. 
Appendix  B  presents  the  Matlab  code. 
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D.  OFDM-BASED  IEEE  802.11 A  STANDARD 


The  rapid  growth  of  the  applications  utilizing  digital  communication  systems  in¬ 
creased  the  need  for  high-speed  data  transmission.  New  multi- carrier  modulation  tech¬ 
niques  are  being  proposed  and  implemented  to  keep  up  with  the  demand  of  higher  data 
rates.  Of  these  multi-carrier  techniques,  OFDM  is  the  method  of  choice  for  high-speed 
communication  due  to  its  many  attractive  features.  This  chapter  attempts  to  justify  the 
choice  of  OFDM  among  other  communication  techniques. 

A.  INTRODUCTION  TO  OFDM 

OFDM  was  introduced  in  the  1950s  but  was  first  implemented  in  the  1960s.  It 
was  originally  developed  from  the  multi- carrier  modulation  techniques  used  in  high  fre¬ 
quency  military  radios.  A  patent  for  OFDM  was  granted  in  the  1970s.  However,  when 
OFDM  was  first  introduced,  it  was  not  very  popular  because  of  the  cost  and  complexity 
of  large  arrays  of  sinusoidal  generators  and  coherence  demodulators  [9]. 

The  actual  widespread  use  of  OFDM  started  after  the  inverse  discrete  Fourier 
transform  (IDFT)  and  discrete  Fourier  transform  (IDFT)  made  the  OFDM  implementa¬ 
tion  possible  without  the  use  of  large  number  of  sinusoidal  generators. 

OFDM  was  accepted  as  a  wireless  local  area  network  (WLAN)  standard  in  Sep¬ 
tember  1999,  but  actually  it  was  not  the  first  IEEE  physical  standard  for  WLANs.  The 
first  standard  was  approved  in  June  1997  and  specified  one  medium  access  control 
(MAC)  and  three  physical  layers  (IEEE  802.11  FHSS,  IEEE  802.11  DSSS  and  IEEE 
802.11  IR).  The  IEEE  802.11  direct  sequence  spread  spectrum  (DSSS)  originally  sup¬ 
ported  both  1  Mbps  and  2  Mbps  of  data  rates  while  the  other  two  supported  1  Mbps  or  2 
Mbps  optionally.  Due  to  the  increasing  demand  for  higher  throughput,  a  high-data-ratc 
DSSS  (IEEE  802. 1  lb)  was  selected  for  standardization  in  July  1998,  and  the  data  rate 
increased  to  1 1  Mbps.  The  IEEE  802. 1  la  and  IEEE  802. 1  lb  standards  were  developed 
simultaneously.  The  federal  communications  commission  (FCC)  released  300  MHz  of 
spectmm  in  the  5.2- GHz  band  in  January  1997  for  WLAN  applications,  and  IEEE 
802.1  la  was  targeted  to  use  this  spectral  band  [10]. 
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OFDM  was  chosen  to  be  the  physical  standard  for  IEEE  802.1  la.  The  time- 
related  parameters  of  OFDM  for  IEEE  802. 1  la  are  listed  in  Table  1. 


Parameter 

Value 

Number  of  data  sub-carriers 

48 

Number  of  pilot  sub-carriers 

4 

Total  number  of  sub -carriers 

52 

Sub-carrier  frequency  spacing 

0.3125  MHz 

IFFT/FFT  period 

3.2  ps  (1/ Af  ) 

Preamble  duration 

16ps 

Signal  Duration  BPSK-OFDM 
symbol 

4  Ms  ( Tgi+Tfft ) 

Guard  Interval  (GI)  duration 

0.8  ps  (W  4) 

Training  symbol  GI  duration 

1.6ps  {Tfft/2) 

Symbol  interval 

4  Ms  {Tgi+Tfft) 

Short  training  sequence  duration 

8ps  (10Tfft/4) 

Long  training  sequence  duration 

8  ms  (Tci+2Tfft) 

Table  1.  OFDM  Timing  Related  Parameters  (From  Ref.  11.) 


The  rate  dependent  parameters  of  IEEE  802.1  la  are  listed  in  Table  2.  As  can  be 
observed  from  Table  2,  different  sub-carrier  modulation  schemes  can  be  used  in  combi¬ 
nation  with  different  convolutional  codes  to  achieve  the  desired  data  rate  for  the  applica¬ 
tion. 


Data  rate 
(Mbps) 

Modulation 

Coding 

Rate(R) 

Coded  bits 
per  sub¬ 
carrier 

(A tbpsc) 

Coded  bits 
per  OFDM 
symbol 

(/Vcsps) 

Data  bits 
per  OFDM 
symbol 

(A  dbps) 

6 

BPSK 

1/2 

1 

48 

24 

9 

BPSK 

3/4 

1 

48 

36 

12 

QPSK 

1/2 

2 

96 

48 

18 

QPSK 

3/4 

2 

96 

72 

24 

16- QAM 

1/2 

4 

192 

96 

36 

16- QAM 

3/4 

4 

192 

144 

48 

64- QAM 

2/3 

6 

288 

192 

54 

64- QAM 

3/4 

6 

288 

216 

Table  2.  OFDM  Rate  Dependent  Parameters  (From  Ref.  11.) 
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Table  3  shows  the  major  parameters  of  the  OFDM  physical  layer.  In  the  simula¬ 
tions  conducted  in  this  thesis,  a  data  rate  of  24  Mbps,  modulation  type  of  QPSK,  and  a 
coding  rate  of  Vi  were  used,  while  the  other  parameters  in  Table  3  were  left  unchanged. 


Information  data  rate 

6,9,12,18,24,36,48  and  54  Mbps 
(6,12  and  24  Mbps  are  mandatory) 

Modulation 

BPSK  OFDM 

QPSK  OFDM 

16- QAM  OFDM 

64- QAM  OFDM 

Error  correcting  code 

K  =  7  convolutional  code 

Coding  rate 

1/2,  2/3,  3/4 

Number  of  sub-carriers 

52 

OFDM  symbol  duration 

4ps 

Guard  interval 

0.8  ps  (Ta) 

Occupied  bandwidth 

16.6  MHz 

Table  3.  Major  Parameters  of  OFDM  Physical  Layer  (From  Ref.  11.) 

OFDM  is  being  used  in  a  number  of  wired  and  wireless  voice  and  data  applica¬ 
tions  due  to  its  flexible  system  architecture.  Some  examples  of  current  OFDM  applica¬ 
tions  are  DAB  (digital  audio  broadcasting),  HDSL  (High- Rate  Digital  Subscriber  Line), 
VHDSL  (Very  High- Rate  Digital  Subscriber  Line),  ADSL  (Asymmetric  DSL),  HDTV- 
Terrestrial  (High  Definition  TV-T),  IEEE  802.11  and  HiperLAN/2,  GSTN  (General 
Switched  Telephone  Network),  Cellular  Radio  and  DVB-T  (digital  video  broadcasting- 
terrestrial)  [9]. 

From  the  range  of  the  current  applications  and  the  research  activity  on  OFDM  and 
considering  the  fact  that  many  wireless  communication  systems  and  network  standards 
are  being  developed  based  on  OFDM,  it  may  be  easily  remarked  that  this  technique  will 
have  a  significant  impact  on  the  future  of  digital  communications  and  wireless  networks. 
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B.  FUNDAMENTALS  OF  OFDM 


1.  Multi-Path  (Delay-Spread  or  Time  Dispersion) 

In  general,  high  data  rate  means  short  symbol  time  compared  to  the  delay  spread 

[j symbol  <  T May )  •  Delay- spread  greatly  affects  the  communication  system  and  the  signal 
might  not  be  recovered  at  the  receiver. 

This  section  addresses  the  effects  of  delay  spread  which  occurs  as  the  surfaces  be¬ 
tween  a  transmitter  and  a  receiver  reflect  a  transmitted  signal  as  in  Figure  1. 


Figure  1.  A  Multi- Path  Scenario 


The  receiver  obtains  the  transmitted  signals  with  random  phase  offsets  and  this 
causes  random  signal  fades  as  reflected  signals  destructively  or  constructively  affect  each 
other  [12],  as  seen  in  Figure  2. 


Transmitted  signal 


r 

Line  of  sight  signal 


Received 

signals  v  Reflected  signal 


Reflected  signal 
V 


Received  symbols  add  up 
at  the  receiver 


Figure  2.  Delayed  Signals  (After  Ref.  13.) 
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When  Tsymhol  < Tdelay  (Bc  <  Bs) ,  as  in  Figure  3,  the  signal  faces  frequency  selective 

fading  and  this  causes  time  dispersion.  The  effect  of  this  is  intersymbol  interference  (ISI), 
where  the  energy  of  one  symbol  leaks  into  another  symbol,  as  can  be  viewed  from  Figure 
4.  As  a  result,  the  bit  error  rate  (BER)  increases,  this  in  turn  degrades  the  performance. 

ISI  is  one  of  the  biggest  problems  of  digital  communication  and  OFDM  deals  with  this 
problem  vciy  effectively. 


Be 


(a)  (b) 

Figure  3.  Representation  of  a  Symbol  in  A  Frequency  Selective  Channel:  (a)  Time 

domain  (b)  Frequency  domain 
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Figure  4.  Illustration  of  ISI  (After  Ref.  14.) 
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A  way  to  deal  with  frequency  selective  fading  is  to  decrease  the  data  rate  and  thus 
change  the  frequency  selective  fading  to  flat  fading.  The  desired  scheme  is  illustrated  in 
Figure  5.  OFDM  systems  mitigate  the  ISI  by  changing  the  frequency  selective  fading 
channel  to  flat  fading  channel  as  discussed  below. 


B: 


Ts-  Syrrtiol  Period 
Td-  Delay  Spread 


_ ^  Time 

Ts 


Figure  5.  (a)  Time  Domain  Representation,  (b)  Frequency  Domain  Representation  of 

a  Symbol  in  Flat  Fading  Channel. 


OFDM  modulates  user  data  onto  tones  by  using  either  phase  shift  keying  (PSK) 
or  quadrature  amplitude  modulation  (QAM).  An  OFDM  system  takes  a  high  data  rate 
stream,  splits  it  into  N  parallel  data  streams  and  transmits  them  simultaneously.  As  can  be 
observed  from  Figure  6,  each  of  these  parallel  data  streams  has  a  rate  of  R/N ,  where  R  is 
the  original  data  rate.  The  data  streams  are  modulated  by  different  carriers  and  combined 
together  by  inverse  fast  Fourier  transform  (IFFT)  to  generate  the  time-domain  signal  to 
be  transmitted  [12]. 
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Figure  6.  OFDM  Splits  a  Data  Stream  into  N  Parallel  Data  Streams  (After  Ref.  15.) 
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By  creating  a  slower  data  stream,  the  symbol  duration  becomes  larger  than  the 
channel’s  impulse  response.  In  this  way,  each  carrier  is  subject  to  flat  fading. 

2.  Orthogonality 

Normally,  multicarrier  systems,  such  as  frequency  division  multiplexing  (FDM), 
have  to  modulate  different  sub -carriers  with  spectrally  separate  symbols  to  avoid  inter 
carrier  interference  (ICI)  at  the  cost  of  a  bandwidth  loss.  However,  in  OFDM,  spectrally 
overlapped  sub -carriers  can  be  used  and  since  they  are  orthogonal,  they  do  not  interfere 
with  each  other.  This  makes  OFDM  a  bandwidth  efficient  modulation  scheme  [12]. 

Orthogonality  of  the  sub-carriers  must  be  ensured  to  avoid  ICI.  OFDM  carriers 
are  orthogonal  over  a  symbol  interval  only  if  they  are  spaced  in  frequency  exactly  at  the 
inverse  of  the  symbol  period  (see  Figures  7  and  8).  IFF  !  implementation  inherently  pro¬ 
vides  this  frequency  spacing  requirement.  It  is  also  necessary  to  add  a  guard  interval  in 
multi- path  channels  to  account  for  the  multi-path  spreading  effects.  Figure  7  shows  a  sub- 
carrier  having  a  rectangular  pulse  shape  and  its  frequency  spectrum.  The  rectangular 
pulse  has  a  width  of  Ts  and  has  a  sine -spectrum  with  nulls  at  multiples  of  l/Ts  . 


Figure  7.  Representation  of  an  OFDM  Sub-Carrier:  (a)  In  Time  Domain  (B)  Fre¬ 

quency  Domain 


Figure  8  shows  the  spectra  of  multiple  carriers.  Although  all  frequency  spectra  are 
overlapping,  they  do  not  overlap  at  the  carrier  frequencies  /,,/2,... ,  and  this  implies  or¬ 
thogonality.  In  the  time  domain  each  carrier  has  an  integer  number  of  periods  within  the 
symbol  time  interval.  Due  to  the  orthogonality,  the  nulls  of  the  side  lobes  coincide  with 
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the  peaks  of  the  main- lobes  of  each  carrier.  By  building  the  receiver  to  sample  at  the  cen¬ 
ter  frequency  of  each  sub -carrier,  it  is  possible  to  obtain  only  the  corresponding  signal’s 
energy,  since  the  side-lobes  of  other  carriers  have  zero  energy  at  that  frequency  [12]. 


Figure  8.  Sub-Carrier  Orthogonality  (After  Ref.  16.) 


Orthogonality  can  be  lost  due  to  two  main  reasons.  The  first  reason  is  the  fre¬ 
quency  errors,  which  are  discussed  in  detail  in  Chapter  in,  and  the  second  is  the  error  in 
symbol  boundary  and  sampling  rate  offset. 

3.  OFDM  Transmitter 

A  block  diagram  of  the  OFDM  transmitter  module  is  presented  in  Figure  9.  Each 
of  the  blocks  is  explained  in  detail  in  the  following  subsections. 


RF 
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Figure  9.  Schematic  Diagram  of  an  OFDM  Transmitter 
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a.  Channel  Coding 

A  sequential  binary  input  data  stream  is  first  encoded  by  the  channel 
coder.  Error  correction  coding  is  important  for  OFDM  systems  used  for  mobile  commu¬ 
nications.  When  channel  coding  is  used  to  improve  its  performance,  OFDM  is  referred  to 
as  coded  OFDM  (COFDM).  The  IEEE  802.11a  standard  uses  convolutional  coding  as 
forward  error  correction  (FEC)  coding,  while  the  convolutional  coder  has  rate  V2,  con¬ 
straint  length  7  and  generator  polynomials  (133,  171)  [11]. 

Table  2  lists  the  rate  and  modulation  scheme  combinations  for  a  desired 
value  of  data  rate.  Higher  data  rates  are  obtained  by  puncturing  the  rate  V2  code.  The 
Convolutional  encoder  used  in  the  IEEE  802.1  la  standard  is  presented  in  Figure  10. 


Input 


► 

Output 

Figure  10.  Convolutional  Encoder  with  Constraint  Length=  7 .  Each  box  represents  a 

shift  register  (After  Ref.  11.) 

b.  Interleaving 

The  encoded  symbols,  out  of  the  convolutional  encoder,  are  next  inter¬ 
leaved  by  a  block  interleaver.  The  block  interleaver  takes  data  from  the  encoder,  per- 
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mutes  it  and  then  outputs  the  rearranged  binary  data.  The  idea  behind  interleaving  is  to 
handle  burst  errors.  The  interleaver  simply  scatters  the  burst  errors  randomly  so  that  the 
receiver  can  easily  correct  the  random  errors. 

Many  interleaving  methods  exist,  but  the  conventional  block  interleaver  is 
widely  used  since  it  is  the  easiest  to  implement.  To  implement  block  interleaving,  a  ma¬ 
trix  that  can  accommodate  all  the  symbols  in  the  transmit  message  matrix  must  first  be 
created.  Next,  the  symbols  from  the  message  matrix  are  read  into  this  intermediate  matrix 
by  rows  and  the  interleaving  effect  is  created  by  reading  these  symbols  out  by  columns. 
The  degree  of  interleaving  is  determined  by  the  dimensions  of  this  intermediate  matrix. 

c.  Symbol  Mapping 

For  the  symbol  mapping  part  of  the  block  diagram,  the  required  data  rate 
for  the  application  will  determine  which  signal  constellation,  i.e.,  which  modulation  tech¬ 
nique  is  to  be  used. 

Following  the  interleaver,  the  rearranged  binary  data  are  mapped  onto  the 
desired  constellation  points.  Figure  11  shows  two  constellation  examples  that  were  listed 
in  Table  2. 


Figure  11.  16- QAM  and  QPSK  Constellations  (From  Ref.  11.) 
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cl.  Differential  Quadrature  Phase  Shift  Keying  (DQPSK)  Modula¬ 
tion 

After  the  symbol  mapping  process,  the  data  can  be  frequency  or  time  dif¬ 
ferential  encoded.  The  Matlab  implementation  details  of  differential  encoding  are  ex¬ 
plained  in  Appendix  A. 


Differential  (noncoherent)  modulation  techniques  decode  the  incoming 
symbols  by  comparing  the  received  symbol  with  the  preceding  symbol.  If  the  transmitted 
waveform  is  given  as  [17] 


4* ) 


cos[w0f+e,.(0] 


i=l,...,M  0  <t<T, 


(0.1) 


then  the  received  signal  is  given  by 


l2j? 

■(t)  =  J—cos[w0t  +  ei(t)+a]  +  n(t) 


i  =  0  <t<T 


(0.2) 


where  a  is  assumed  to  be  a  random  variable  uniformly  distributed  between  0  and  2 %  and 
n  ( t )  is  an  additive  white  Gaussian  noise  (AWGN)  process. 


For  noncoherent  detection,  matched  filters  cannot  be  used  because  the 
matched  filter  output  becomes  a  function  of  the  unknown  phase  a  .  However,  assuming 
that  a  changes  slowly  relative  to  the  two  symbol  period  times,  the  phase  difference  of 
two  respective  signals  becomes  independent  of  a  and  given  by  [17] 

[e1(r2)+a]-[e;(r)+a]=et(r2)-eJ(r1)=(|l,(7-2),  (o.3> 

which  suggests  that  the  phase  of  the  previous  signaling  interval  can  be  used  as  a  phase 
reference  to  demodulate  the  current  signal.  Of  course,  this  requires  the  use  of  differential 
encoding  in  the  transmitter. 

Even  though  DPSK  performs  less  efficiently  than  PSK,  it  was  chosen  due 
to  the  ease  of  implementation. 

e.  IF  FT 

Hie  interleaved  and  differentially  encoded  complex  frequency  array  is 
next  modulated  onto  a  desired  number  of  sub -carriers  (52  sub -carriers  are  used  in  the 
TREE  802.1  la  standard)  by  implementing  an  IFFT. 
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The  OFDM  baseband  sub -carrier  is 

<t>*  (0  =  ej2%fk>  (0.4) 

where  fk  is  the  k ,h  sub -carrier  frequency.  An  OFDM  symbol  consists  of  iV  modulated 
sub-carriers.  The  OFDM  signal  not  including  a  cyclic  prefix  is  given  by  [3] 

0  <t<  NT  (0.5) 

k=  0 

where  z,k is  the  k'h  complex  data  symbol  and  NT  is  the  OFDM  symbol  duration.  The 
sub-carriers  in  Equations  2.4  and  2.5  have  frequencies 


in  the  sense  that 


(0.6) 


,  NT 

™  =  (0.7) 

iVi  0 

ensures  orthogonahty.  If  the  signal  s  (t )  is  sampled  with  a  sampling  period  of  T,  the  fol¬ 
lowing  is  obtained: 

s{n)  =  s{t  =  nT)  =  \D¥T{zk}  =  -^rYJzkei27lknlN  «  =  0,1...A-1.  (0.8) 

This  last  equation  is  IDFTjz,  ]  and  was  proposed  by  [18].  As  can 

be  seen  from  Equation  2.8,  a  baseband  OFDM  transmission  symbol  is  an  iV- point  com¬ 
plex  modulation  sequence.  It  is  composed  of  N  complex  sinusoids,  which  are  modulated 
with  z(k) . 

f.  Guard  Interval  Addition 

Hie  addition  of  the  Guard  interval,  TCI ,  is  the  next  step  in  generating 
OFDM  symbols.  The  guard  interval  is  added  to  the  start  of  OFDM  symbols  at  the  trans¬ 
mitter  to  eliminate  ISI  and  ICI,  and  it  is  deleted  before  the  EFT  process  at  the  receiver. 


Hie  interval  TC1  is  an  important  parameter  in  the  IEEE  802.1  la  standard 
because  it  determines  the  choice  of  other  parameters.  It  is  normally  desired  that  TC1  be 
larger  than  the  expected  multi-path  delay  spread  to  combat  ISI.  However,  the  larger  Tcl , 


14 


the  smaller  the  effective  symbol  duration,  thereby  resulting  in  reduced  throughput.  In  the 
IEEE  802.1  la  standard,  the  total  OFDM  symbol  duration  is  chosen  to  be  4  ps  ,  and  TCI  is 
0.8  ps  .  Therefore,  the  effective  symbol  duration  is  3.2  ps  ,  and  the  inverse  of  the  effec¬ 
tive  symbol  duration  provides  the  sub-carrier  spacing  of  0.3125  MHz.  These  time  related 
parameters  of  OFDM  are  listed  in  Table  1. 

The  expression  for  a  low- pass  OFDM  signal  is  given  by  Equation  2.8.  The 
transmitted  time  domain  signal  is  then  obtained  by  including  the  guard  interval.  The  ana¬ 
log  signal  after  the  digital  to  analog  converter  (DAC)  in  the  block  diagram  is  expressed  as 

(0 

A  k=0 

where  A/  is  sub- carrier  spacing.  The  guard  interval  Tqi  is  created  by  appending  the  IFFT 
sequence  with  a  circular  extension  (called  cyclic  prefix)  in  order  to  maintain  orthogonal¬ 
ity.  For  example,  the  cyclic  prefix  for  the  transmission  sequence  j  slp  (/?)}  consists  of 

the  last  m  samples  { slp ( N-m),s Ip ( N  - m  - 1 ) , . . .  sjp(N  -1)} ,  which  are  added  to  the  begin¬ 
ning  of  the  sequencej  ^(0),  s,  (l),....?,  (iV-1)  j .  The  new  sequence  now  consists  of  N+m 

samples.  The  benefit  of  the  guard  interval  can  be  observed  in  Figure  12.  Figure  12b  illus¬ 
trates  that  the  decaying  transient  of  symbol  1  in  the  guard  interval  can  be  allowed  since  it 
does  not  affect  symbol  2. 
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(b)  OFDM  Symbols  with  Guard  Interval 
Figure  12.  Effects  of  Guard  Interval  (After  Ref.  16.) 


The  disadvantage  of  the  guard  interval  is  that  it  causes  a  loss  of  bandwidth 
efficiency  since  data  are  not  transmitted  during  the  guard  interval.  On  the  other  hand,  by 
adding  the  guard  interval,  ICI  can  be  avoided  since  orthogonality  is  ensured,  and  ISI  is 
avoided  since  adjacent  symbols  do  not  interfere  with  each  other. 
g.  Number  of  Sub-carriers 

The  number  of  sub -carriers  is  another  important  parameter  in  the  design  of 
OFDM  systems  along  with  the  system  bandwidth  B  ~  l/T  ,  the  sub -carrier  bandwidth 

1/ NT  ,  and  the  guard  interval  TCI . 

Initially,  the  guard  interval  TC1  is  chosen  to  be  larger  than  delay  spread  x  , 
and  the  OFDM  symbol  length  should  be  larger  than  TC1 .  Therefore  NT  »  x  ,  where  N  is 

the  number  of  sub -carriers  which  yields  N  »  Bz  .  However,  a  large  N  means  smaller  fre¬ 
quency  spacing  between  OFDM  carriers,  which  in  turn,  makes  the  system  vulnerable  to 
Doppler  spread  and  ICI.  This  last  restriction  forces  larger  frequency  spacing  than  the 
Doppler  spread  (B/N  »  fd ) .  As  a  result,  the  number  of  sub-carriers  can  be  chosen  in  the 
interval  given  by  [3] 
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The  number  of  sub-carriers  is  fixed  at  52  in  the  TERR  802.1  la  standard.  Pour  of 
these  sub-carriers  are  used  as  a  reference  to  rrunirnize  frequency  and  phase  shifts  of  the 
signal  during  transmission. 

h.  Symbol  Pulse  Shaping  ( Windowing ) 

The  time  dispersion  can  be  combated  by  the  addition  of  guard  interval. 
However,  it  is  also  necessary  to  combat  the  slow  decay  of  OPDM  symbols  in  the  fre¬ 
quency  domain.  The  effects  of  slow  decay  into  the  adjacent  bands  can  be  avoided  by  nar¬ 
rowing  the  spectrum. 

The  method  of  choice  to  narrow  the  spectra  is  to  use  pulse  shaping  of  the 
OPDM  symbol.  Pulse  shaping  is  accomplished  by  using  either  a  time  window  or  a  filter 
[3].  Windowing  helps  to  smooth  the  transition  between  symbols  and  narrows  the  signal 
spectrum. 

i.  RF  Modulation 

After  the  pulse  shaping  process,  the  OFDM  symbols  are  upconverted  to  a 
radio  frequency  carrier,  amplified  and  transmitted  through  the  antenna.  The  TPPP 
802.1  la  standard  specifies  the  5-Ghz  band  for  the  carrier  frequency.  In  this  thesis,  the 
communication  system  was  simulated  at  the  baseband. 

4.  OFDM  Receiver 

A  block  diagram  of  the  OFDM  Receiver  module  is  presented  in  Figure  13. 


RF 

demodulation 


Figure  13.  Schematic  Diagram  of  an  OFDM  Receiver 
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The  processes  performed  in  the  transmitter  are  simply  reversed  in  the  receiver.  In 
addition,  the  receiver  must  include  algorithms  for  synchronization  and  channel  estima¬ 
tion.  In  particular  the  synchronization  block  performs  operations,  such  as  packet  detec¬ 
tion,  timing  estimation,  frequency  offset  estimation  and  correction. 

a.  Removing  Guard  Interval  and  FFT  Processing 

Since  only  the  baseband  processing  is  being  used  in  this  thesis,  the  down- 
conversion  process  in  the  receiver  is  not  needed.  The  first  step  is  to  remove  the  guard  in¬ 
terval  to  obtain  the  information  portion  of  the  symbol  for  further  processing.  Next,  the 
time  domain  samples  are  transformed  into  the  frequency  domain  by  the  FFT  process. 

This  also  makes  it  possible  to  recover  the  OFDM  frequency  tones. 

b.  Decoding  and  Deinterleaving 

The  next  step  in  the  receiver  is  the  time  or  frequency  differential  decoding. 
Following  the  differential  decoding,  the  inverse  mapping  of  each  received  complex 
modulation  value  into  a  corresponding  N-  ary  symbol  is  accomplished.  The  message  is 
next  deinterleaved  according  to  the  particular  interleaving  configured  in  the  transmitter. 

c.  Received  message  and  Viterbi  Decoding 

After  being  exposed  to  the  effects  of  different  channel  models,  some  errors 
may  occur  in  the  transmitted  message.  The  received  bit  stream  is  determined  after  the 
Viterbi  decoder. 

The  Viterbi  decoders  are  used  to  decode  bit  streams  encoded  by  convolu¬ 
tional  encoders.  After  receiving  the  bit  stream  from  the  channel,  the  Viterbi  algorithm 
searches  the  trellis  to  find  the  path  that  can  create  the  bit  stream  closest  to  the  received 
one.  The  measure  of  closeness  is  different  for  hard  and  soft  decision  decoding.  For  the 
hard  decision  decoding,  the  Viterbi  algorithm  tries  to  find  the  bit  stream  path  that  has  the 
minimum  Hamming  distance  from  the  received  sequence.  For  soft  decision  decoding,  the 
algorithm  tries  to  find  the  bit  stream  path  that  has  the  minimum  Euclidian  distance  from 
the  received  sequence.  The  details  of  the  Viterbi  decoding  algorithm  can  be  found  in  [17]. 
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C.  SUMMARY 

An  introduction  to  OFDM  was  presented  in  this  chapter.  The  transmitter  and  re¬ 
ceiver  block  diagrams  have  been  discussed  in  detail  and  the  chapter  concluded  by  listing 
the  advantages  and  disadvantages  of  OFDM  and  mentioning  the  potential  problematic 
areas  in  OFDM  implementation. 

The  advantages  of  OFDM  can  be  summarized  as  follows.  OFDM  changes  fre¬ 
quency  selective  fading  to  flat  fading  and  offers  high  bit  rate  transmission.  Its  flexible 
architecture  allows  the  use  of  combinations  of  various  coding  schemes  (rates)  and  modu¬ 
lation  schemes  (signal  constellations)  according  to  channel  condition.  OFDM  eliminates 
ISI,  is  bandwidth  efficient,  and  supports  several  multiple  access  schemes  (TDMA, 

FDMA,  MC-CDMA,  etc.)  [15]. 

On  the  other  hand  OFDM  is  very  sensitive  to  frequency  errors,  has  a  large  peak  to 
average  ratio  (PAR),  and  is  sensitive  to  channel  fading  [15].  The  effect  of  PAR  is  ana¬ 
lyzed  in  [19]. 

OFDM  has  some  potential  problematic  areas  that  should  be  considered  during 
implementation.  The  first  is  synchronization,  including  the  time  and  frequency  synchro¬ 
nization.  The  second  is  the  non-constant  power  envelope,  which  requires  the  use  of  linear 
amplifiers.  The  last  problematic  area  is  the  requirement  of  channel  estimation  since  the 
channel  is  time- variant  [13]. 

The  next  chapter  presents  the  steps  to  create  an  OFDM  signal.  OFDM  systems’ 
unique  requirements,  such  as  IFFT/FF1  algorithms  and  the  cyclic  prefix,  are  also  cov¬ 
ered.  Following  that,  the  effects  of  frequency  errors  are  analyzed. 
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ID.  FREQUENCY  ERRORS  IN  OFDM  AND  THEIR  EFFECTS 


Before  an  OFDM  receiver  demodulates  the  sub-carriers,  it  performs  two  synchro¬ 
nization  tasks.  First,  it  determines  the  symbol  boundaries  and  the  optimal  timing  instants 
in  order  to  minimize  both  the  ICI  and  the  ISI.  Second,  it  tries  to  estimate  and  correct  the 
frequency  errors  [20].  OFDM’s  sensitivity  to  these  frequency  errors  is  one  of  its  main 
drawbacks.  This  chapter  discusses  the  effects  of  these  frequency  errors  on  OFDM  sys¬ 
tems. 

A.  FREQUENCY  ERRORS 

The  orthogonality  of  the  sub -carriers  can  be  ensured  only  if  the  receiver  and  the 
transmitter  have  the  same  reference  frequency.  Any  deviation  from  this  reference  fre¬ 
quency  may  cause  ICI  and  loss  of  orthogonality.  Another  frequency  error  factor  is  phase 
noise,  which  is  caused  by  random  jitter  of  the  phase  of  the  steady  sinusoidal  waveform 
generated  by  the  oscillators  [20]. 

Typically  frequency  errors  are  generated  by  the  fact  that  the  oscillators  in  the 
modulator  and  demodulator  do  not  have  exactly  the  same  frequency. 

For  single- carrier  systems,  the  effect  of  phase  noise  and  frequency  offset  appear 
only  as  degradation  in  the  received  SNR,  rather  than  ISI  or  ICI.  Nevertheless,  many  effi¬ 
cient  techniques  to  minimize  the  effects  of  this  drawback  have  been  proposed  in  the  lit¬ 
erature  [20,  21]. 

Other  reasons  for  frequency  errors  include  Doppler  shift  caused  by  the  relative 
movement  between  the  receiver  and  the  transmitter,  and  phase  noise  introduced  by  non¬ 
linear  channels. 

Figure  15  shows  the  front  end  of  an  OFDM  receiver  where  most  of  the  frequency 
errors  occur,  i.e.,  the  local  oscillators  and  the  sample  clock  at  the  analog  to  digital  (A/D) 
converter  [3].  The  A/D  converter  causes  errors  when  the  receiver  does  not  have  the  same 
sample  clock  frequency  as  at  the  transmitter. 
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Figure  14.  OFDM  Receiver  Front  End  (After  Ref.  3.) 


The  frequency  offset  and  the  phase  noise  cause  a  phase  rotation  of  the  received 
symbols.  Coherent  OFDM  systems  need  a  phase  tracking  device  to  obtain  the  phase  of 
the  incoming  symbols  for  correct  demodulation  [22].  Typically,  three  types  of  algorithms 
are  used  to  estimate  the  frequency  offsets,  i.e.,  to  track  the  phase,  in  coherent  OFDM  sys¬ 
tems: 

•  data- aided  algorithms  that  use  special  training  information  (pilot  tones) 
within  the  transmitted  signal, 

•  non- data- aided  algorithms  that  analyze  the  received  signal  in  the  fre¬ 
quency  domain,  and 

•  algorithms  that  use  the  cyclic  prefix  of  the  OFDM  signals  [21]. 

For  wireless  applications,  the  first  method  is  the  method  of  choice.  The  imple¬ 
mentation  details  of  these  coherent  OFDM  detection  algorithms  can  be  found  in  [21].  It 
should  be  noted  that  in  many  applications  as  well  as  in  this  thesis,  the  differential  detec¬ 
tion  technique  is  used  since  it  does  not  require  channel  estimation  and  the  use  of  pilot 
tones,  even  though  it  causes  a  loss  of  SNR. 

In  the  following  sections,  the  effects  of  each  of  the  aforementioned  reasons  that 
cause  frequency  errors  are  analyzed.  A  simple  representation  of  an  OFDM  signal  is  cre¬ 
ated  for  the  ease  of  analysis. 
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B.  CREATING  A  SIMPLE  REPRESENTATION  OF  AN  OFDM  SIGNAL  TO 
SIMPLIFY  ANALYSIS 

A  baseband  OFDM  signal  can  be  represented  by  [7] 

N-\ 

b(t)  =  Y;Acos(,X>lt  +  §,)  (0.11) 

;= 1 

where  Ai  is  the  amplitude,  co,  =  2n  f]  is  the  angular  frequency,  A  is  the  phase  of  the  ith 
sub -carrier,  and  N  is  the  number  of  sub -carriers. 


According  to  the  modulation  technique  to  be  used,  either  A  or  (|)  is  determined  by 
the  data.  As  discussed  in  Chapter  n,  the  baseband  OFDM  signal  b(t)  is  modulated  next, 
onto  a  RF  carrier  with  frequency  fc : 


s(t )  =  2b(t)cos(Oct 

N- 1 

=  2^  A  cOS(CQ,.t  +  A)coS  (Dct 

1=0 

N-\ 

=  £  A  {cos  [K  +©,■  )t+A]  +c°s  [(©„  -co, )  t  -  A  ]} 

1=0 


(0.12) 


where  ot  =  2n  /'  ,  and  we  assume  the  phase  of  the  carrier  to  be  zero  for  simplicity.  Since 
a  single  side  band  transmission  is  enough  to  carry  the  information  in  A  or  A  ,  it  is  as¬ 
sumed  that  the  upper  sideband  is  used,  and  therefore  the  transmitted  signal  can  be  repre¬ 
sented  as  [7] 

N- 1 

s(t)  =  Y,  A  cos  [(co,  +  <o,  )t +<!>,']•  (0.13) 

i=0 

C.  THE  EFFECT  OF  DOPPLER  FREQUENCY  SHIFT,  PHASE  NOISE  AND 
OSCILLATOR  FREQUENCY  OFFSET  ON  THE  OFDM  SYSTEMS 

In  this  section  the  theoretical  analysis  of  the  effects  of  frequency  errors  is  pre¬ 
sented.  The  maximum  Doppler  shift  occurs  when  the  two  mobile  nodes  move  toward 
each  other,  given  by  [23] 
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c 

where  v  is  the  relative  speed  of  the  two  nodes,  fc  is  the  carrier  frequency  and  c  is  the 


speed  of  light(3xl08m/s) . 

1.  Effect  of  Doppler  Shift  on  the  Carrier  Frequencies,  Sub-carriers,  En¬ 
velope  and  Symbol  Timing 

An  OFDM  signal  consists  of  numerous  sub-carriers  with  different  frequencies. 

The  amount  of  Doppler  shift  affecting  the  i'h  sub-carrier  is  given  by  [7] 

(Jc±fi)  DOPPler  m  >(l  +  0(/c±/|)  (°-15) 

where  c,  is  the  percentage  of  the  change  in  frequency  and  is  determined  by 

£,  =  — =  — cos0.  (0.16) 

/  c 

Tlie  right-hand  side  of  Equation  3.5  can  be  written  as 

(i+5)(/c±/,)=(i+5)/c±(i+5)/,  (0-17) 

which  demonstrates  that  the  Doppler  frequency  shift  affects  the  carrier  frequency  and  the 
sub -carrier  frequencies  by  the  same  percentage  c,  [7].  The  Doppler  shift  of  the  carrier 

frequency  can  be  calculated  as 

fdc  =  — cose  (0.18) 

c 

and  the  Doppler  shift  of  the  sub -carrier  frequencies  as 

fa=~  cosO.  (0.19) 

c 

By  using  Equation  3.3  again,  the  transmitted  OFDM  signal  with  Doppler  shift  can 
be  written  as 


N- 1 


s(t)  =  £a,  cos[(1+S)( cq  +C0,  )t+f.] 


/= o 

N- 1 


(0.20) 


=  E { A  cos  [( 1  ]  cos [  ( 1  +  £ )  0)/]  -  A,,  sin  [( 1- +  £ ) co/  +4>,  ]  sin  [(1- +  £ ) co/]}. 

i= 0 
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In  Equation  3.10,  A.  cos  [( 1  +  q  )cqf  +  (j);]  can  be  thought  of  as  the  envelope  of  the 

carrier,  cos  [( 1  +  £, )  co  t J  ,  which  helps  to  demonstrate  that  the  Doppler  shift  affects  the 

envelope  and  the  carrier  frequency  by  the  same  percentage  [7].  The  Doppler  shift  also 
affects  the  symbol  rate  and  the  time  synchronization. 

2.  Phase  Noise,  Oscillator  Frequency  Offset  and  Their  Effects 
The  carrier  frequency  offset  has  two  destructive  effects  on  OFDM  systems.  The 
first  is  the  reduction  of  the  amplitude  of  the  desired  sub -carrier  and  the  second  is  the  in¬ 
troduction  of  ICI  [21].  The  reduction  of  the  amplitude  happens  because  the  desired  sub- 
carrier  is  not  sampled  at  the  peak  of  the  sine  function  of  the  DFT  since  the  sine  functions 
are  shifted  (see  Figure  15).  Also,  the  ICI  is  caused  since  orthogonality  is  lost  between  the 
neighboring  sub -carriers  [24]. 


Frequency 


Moose  [5]  and  Pollet  et  al.  [4]  analytically  evaluated  the  effects  of  the  earner  fre¬ 
quency  offset  and  carrier  phase  noise  on  the  SNR  degradation  for  an  AWGN  channel. 
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The  results  derived  in  these  papers  are  used  for  many  of  the  studies  on  frequency  errors. 

As  a  result,  we  have  used  these  results  as  benchmarks  for  the  simulation  results  in  this 
thesis. 

Moose  [5]  models  the  channel  by  a  complex  transfer  function  Hk  for  the  k'b  sub- 
carrier  and  assumes  that  a  frequency  offset  of  £  exists  for  all  sub-carriers.  The  frequency 
offset,  £  ,  is  usually  expressed  relative  to  the  data  symbol  rate,  and  it  includes  both  the 
Doppler  shift  and  the  local  oscillator  frequency  offset  effects.  As  mentioned  previously, 
each  sub-carrier  is  affected  by  the  Doppler  shift  according  to  its  frequency  value.  Taking 
all  the  causes  into  consideration,  the  total  offset  for  each  sub -carrier  after  down  conver¬ 
sion  can  be  expressed  as  A =  fdc+  fdi  +  A fL0  where  fdc  is  the  carrier  Doppler  shift,  fdi  is 

the  i'h  sub-carrier  Doppler  shift  and  A fLO  is  the  local  oscillator  frequency  offset  [7]. 


As  can  be  seen  from  Equations  3.8  and  3.9,  fdc  »  fdi  since  fc  »  f] .  Thus, 

Afi  =  fdc  +  A fL0  =  A/  .  This  result  demonstrates  that  the  frequency  offset  is  independent  of 
the  sub-carriers,  which  in  turn  states  that  the  relative  offset  £  =  A f  /Rs  is  also  independ¬ 
ent  of  sub-carriers  [7]. 

The  SNR  calculated  at  the  output  of  the  DFT  of  the  receiver  is  given  by  [5] 


SNR  > 


N„ 


sin  TIE 
7t£ 


1 


J  1  +  0.5947 


(sin7t£  )' 


|e|<0.5 


where  Ec/Na  is  the  OFDM  carrier  energy  to  the  AWGN  spectral  density.  An  upper 
bound  on  the  degradation  can  be  obtained  from  Equation  3.1 1  as  follows: 


A*  <101og 
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f  E  ,  3 

1+  0.5947  -^sin27tA/ 
Np 

sin  c2Af 


(0.21) 


(0.22) 


V  ) 

where  the  factor  0.5947  is  derived  by  taking  the  lower  bound  of  the  summation  of  all  in¬ 
terfering  sub -carriers  [24]. 
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Similar  results  to  Equations  3.11  and  3.12  are  found  in  [4]  in  which  the  channel  is 
modeled  by  a  time-varying  phase  0  (t) ,  which  is  a  result  of  either  the  phase  noise  of  the 
carriers  or  a  carrier  offset  between  the  receiver  and  transmitter.  A  distinction  is  made  be¬ 
tween  phase  noise  and  frequency  offset. 

a.  Phase  Noise 

For  the  phase  noise  case,  0  (t)  is  assumed  to  be  a  Wiener  process  with 
£"{0  (/)}  =  0  and£j(0  (to  +  t)-Q  (to  ))"j  =  47tfl  |f|,  where  (Hz)  is  the  one-sided  3-dB 
linewidth  of  the  Lorentzian  power  density  spectrum  of  the  free- running  carrier  generator 

[4]. 


Phase  noise  has  two  main  effects.  First,  it  causes  a  random  phase  variation 
common  to  all  sub -carriers.  The  effects  of  this  common  phase  error  are  minimized  by 
employing  phase  tracking  techniques  or  differential  decoding.  Second,  it  introduces  ICI. 
Based  on  the  model  defined  in  [4],  the  degradation  D  in  SNR,  i.e.,  the  required  increase 
in  SNR  to  compensate  for  the  phase  noise  is 
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Since  R=  N/T  =  NRS ,  where  N  is  the  total  number  of  sub -carriers  and  Rs  is  the  sub- 
carrier  symbol  rate,  Equation  3.13  can  be  rewritten  as 


(0.23) 
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The  plot  of  SNR  degradation  in  dB  as  a  function  of  the  normalized  band¬ 
width  ( /  R  )  is  shown  in  Figure  18  for  QPSK.  This  figure  is  plotted  for 

EjN0  =  10.5  dB  corresponding  to  a  BER  of  1 0  6  for  uncoded  QPSK  [20]. 


(0.24) 
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-3  dB  linewidth/subcarrier  spacing  ^g-3 

Figure  16.  SNR  Degradation  Versus  -  3-dB  Bandwidth  of  the  Phase  Noise  Spectrum 

for  QPSK. 

By  observing  Figure  16,  it  can  be  concluded  that  for  a  negligible  SNR 
degradation  of  about  0.1  dB,  the  -  3-dB  phase  noise  bandwidth  should  be  approximately 
0.1  percent  of  the  sub -carrier  spacing  for  QPSK.  However,  this  value  changes  depending 
on  the  modulation. 

b.  Frequency  Offset 

For  the  frequency  offset  case,  0  it)  is  given  by 

0  (t)  =  2nAFt+Q0  (0.25) 

where  A F  is  the  carrier  offset.  If  there  is  a  frequency  offset,  then  the  number  of  cycles  in 
the  FFT  interval  is  not  an  integer  and  this  causes  ICI  after  the  FFT.  The  sub -carriers  in 
the  middle  of  the  OFDM  spectrum  are  more  exposed  to  ICI  than  the  sub -carriers  on  the 
edges  of  the  spectrum  because  they  have  interfering  sub-carriers  on  their  both  sides  [20]. 
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Based  on  the  model  developed  in  [4],  the  degradation  D  in  SNR  for  fre¬ 
quency  offset  is  given  by 

10 


DAa  =■ 


dB 


3-lnl0 
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By  using  the  relation R  =  N/T  =  NRS ,  Equation  3.16  can  be  rewritten  as 


Aib  —  ' 
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where  A F/Rs  is  the  relative  frequency  offset.  The  plot  of  SNR  degradation  in  dB  as  a 
function  of  the  relative  frequency  offset  (A F/Rs )  is  shown  in  Figure  17  for  QPSK.  This 
figure  is  plotted  for  Es/N0  =  10.5  dB  corresponding  to  a  BER  of  10  6  for  uncoded  QPSK 
[20], 


(0.27) 


Figure  17.  SNR  Degradation  Versus  the  Normalized  Frequency  Offset  for  QPSK. 
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From  Figure  17,  it  can  be  seen  that  for  a  negligible  SNR  degradation  of 
about  0.1  dB,  the  tolerable  frequency  offset  should  be  approximately  2.5%  of  the  sub- 
carrier  spacing  for  QPSK.  However,  this  value  changes  depending  on  the  modulation. 

D.  SUMMARY 

This  chapter  discussed  the  frequency  errors  and  their  causes.  It  was  demonstrated 
that  SNR  degradation  due  to  the  Doppler  shift,  most  of  the  time,  is  negligible,  and  this 
left  frequency  offset  and  phase  noise  of  local  oscillators  as  the  main  sources  of  degrada¬ 
tion.  It  was  also  shown  that  this  degradation  varies  strongly  according  to  the  modulation 
used.  In  [21],  it  was  illustrated  that  64- QAM  could  not  tolerate  more  than  1%  carrier  fre¬ 
quency  error  for  a  SNR  degradation  of  0.5  dB,  whereas  QPSK  modulation  could  tolerate 
5%  carrier  frequency  error  for  the  same  SNR  degradation.  It  can  also  be  stated  that  close 
spacing  of  carriers  in  frequency  made  the  tolerable  frequency  offset  a  small  percentage  of 
the  channel  bandwidth.  In  [5],  it  is  illustrated  by  simulation  that  the  carrier  frequency  off¬ 
set  is  limited  to  4%  or  less  of  the  intercanier  spacing  to  maintain  signal- to- interference 
ratios  of  20  dB  or  greater  for  the  OFDM  carriers. 

In  the  next  chapter,  this  thesis  will  attempt  to  simulate  an  OFDM  system  under 
frequency  offset  and  phase  noise  conditions  to  observe  the  expected  performance  degra¬ 
dation  discussed  in  this  chapter. 
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IV.  OFDM  SYSTEM  SIMULATIONS  IN  MATLAB 


The  previous  chapters  covered  the  principles  of  OFDM  and  discussed  the  fre¬ 
quency  errors  that  might  occur  in  a  communication  system  as  well  as  the  effects  of  these 
frequency  errors.  This  chapter  presents  simulation  of  an  OFDM  communication  system, 
operating  under  different  channel  conditions.  The  results  of  simulation  are  compared  to 
the  theoretical  results  discussed  in  Chapter  IH.  The  simulation  steps  performed  are  pre¬ 
sented  in  Figure  18. 


Figure  18.  Matlab  Simulations  Scheme 

A.  SELECTING  SYSTEM  PARAMETERS 

Hie  simulations  attempt  to  emulate  the  bit  rate  and  performance  requirements  of 
the  IEEE  802.1  la  standard  by  using  the  Matlab  software.  An  operational  bit  rate  of  24 
Mbps  was  chosen  and,  according  to  the  standard,  the  occupied  bandwidth  was  taken  as 
16.6  MHz.  The  constellation  to  be  used  was  decided  based  upon  these  parameters.  First, 
the  number  of  bits/Hz  needed  to  be  transmitted  to  meet  the  requirements  is  calculated  as 

bit.  =  24Mbps  =  j  44bit/Hz  „  2bit/Hz. 

Hz  16.6MHz 
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QPSK  signaling  was  chosen  for  this  simulation,  and  the  actual  bit  rate  was 

2x16.6  =  33 .2Mbps .  Even  though  this  seems  much  more  than  the  required  bit  rate,  the 

addition  of  the  guard  interval  and  EEC  reduces  the  bit  rate  to  the  desired  level  [25]. 

The  required  symbol  interval  in  the  standard  was  given  as  4  ps  ,  and  the  sub- 
canier  spacing  was  given  in  Table  1  as  0.3 125MHz  .  The  number  of  sub-carriers  used  in 

the  standard  is  calculated  using  this  sub -carrier  spacing  and  occupied  bandwidth  as  fol¬ 
lows: 

Ar  occupied  BW  16.6MHz 

N- - - - = - =  53. 

frequencyspacing  0.3125MHz 

The  standard  actually  uses  52  sub -carriers  with  4  of  these  reserved  for  pilot  sig¬ 
nals  in  coherent  detection  systems  to  make  the  system  robust  against  frequency  errors. 
Since  differential  encoding  and  decoding  were  used  in  this  work,  all  52  sub -carriers  could 
be  used  to  carry  information  symbols. 

A  16- point  guard  interval  was  used  during  the  simulations  to  represent  an  800- ns 
period,  and  a  64- point  FFT  was  used  to  represent  the  3.2  ps  information  signal  length, 
i.e.,  a  total  of  80  time  samples  representing  the  standard  4  ps  OFDM  symbol. 

B.  SIMULATION  METHODOLOGY 

After  modifying  the  code  from  [25,  26]  and  adding  required  sub-blocks  to  imple¬ 
ment  the  OFDM  system  for  this  thesis,  a  number  of  system  simulations  were  performed 
using  different  channel  models. 

The  simulation  process  was  performed  in  five  steps.  The  first  step  was  to  verify  if 
the  Matlab  code  was  properly  functioning.  After  verification  of  the  code,  a  channel  with 
only  AWGN  was  simulated,  and  the  system’s  performance  in  this  channel  was  investi¬ 
gated.  The  third  step  was  to  investigate  the  performance  of  four  mobile  channels.  During 
this  step,  how  different  mobile  channels  affect  the  system  was  demonstrated  and,  then, 
the  most  severe  of  them  was  chosen  as  the  mobile  channel.  The  first  three  simulation 
steps  were  conducted  to  ascertain  the  system’s  response  to  different  channel  models 
without  introducing  any  frequency  offset  and  phase  noise.  The  main  goals  of  the  simula¬ 
tions  were  achieved  in  the  fourth  and  the  fifth  steps.  The  fourth  step  included  the  fre- 
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quency  offset,  which  was  studied  for  three  different  channels.  In  the  fifth  step,  the  phase 
noise  was  introduced  and  the  performance  measured  for  three  different  channels.  Table  4 
summarizes  all  of  the  simulation  steps. 


Step  No. 

Simulation 

Channel  Description 

1. 

Code  check 

Channel  without  AWGN  and 

multi- path  effects  (Ideal  channel). 

2. 

Performance  in  AWGN 

AWGN  channel. 

3. 

Performance  in  Mobile 

channel 

Multi- path+ AWGN  channel. 

4. 

Performance  with  Fre¬ 
quency  offset 

a.  Ideal  channel. 

b.  AWGN  channel 

c.  Mobile  channel 

5. 

Performance  with  Phase 

noise 

a.  Ideal  channel. 

b.  AWGN  channel 

c.  Mobile  channel 

Table  4.  Simulation  Steps 


This  step-by-step  simulation  methodology,  beginning  with  the  individual  channel 
simulations  and  advancing  to  the  combined  channel  simulations  helped  to  evaluate  each 
channel  model’s  output  individually.  This  approach  made  it  possible  to  detect  the  errors 
caused  by  incorrectly  designed  sub -blocks  at  an  early  stage  and  forced  a  redesign  of  the 
incorrectly  functioning  blocks. 

The  implementation  details  and  results  of  each  simulation  step  are  explained  be¬ 
low. 
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1.  Code  Check  (Channel  Model  0) 

The  first  step  in  the  simulation  process  was  to  check  whether  or  not  the  system 
model  functioned  satisfactorily  in  a  noise- free,  ideal  channel.  Due  to  the  lack  of  any 
channel  distortions  affecting  the  signal,  the  transmitted  and  received  messages  should  be 
exactly  the  same  in  the  case  of  a  correctly  implemented  overall  system.  After  many  simu¬ 
lations,  with  different  parameters  each  time,  it  was  determined  that  the  system  was  func¬ 
tioning  correctly  since  it  was  possible  to  receive  the  transmitted  signals. 

Table  5  lists  the  parameters  of  a  simulation  run  for  QPSK  realization.  The  trans¬ 
mitted  and  received  QPSK  signal  constellations  under  noise- free  channel  conditions  are 
shown  in  Figure  19,  which  verified  a  correctly  functioning  system  since  we  were  able  to 
receive  exactly  what  we  had  transmitted. 


#  of  car¬ 
riers 

#  of  sym¬ 
bols 

Channel  model 

seed 

Interleaver  ma¬ 
trix 

48 

300 

0 

22 

[28  12] 

Table  5.  Code  Check  Simulation  Parameters 
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Figure  19.  (a)  Transmitted  and  (b)  Received  QPSK  Constellations  After  Differential 

Decoding,  from  a  Channel  Model  0  Simulation 
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2.  OFDM  Performance  in  an  AWGN  Channel  (Channel  Model  1) 

In  this  step,  AWGN  channel  simulations  were  conducted,  and  the  system’s  per¬ 
formance  results  were  compared  with  the  reported  values  in  other  related  studies  [25,  26 
and  27]. 


The  theoretical  symbol  error  rate  Pb  for  BPSK  is  given  by 


Ph  =  Q 


(0.30) 


where  [17] 


Q{x)=^exp 
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For  M-ary  modulation,  Ph  can  be  approximated  by  assuming  that  each  symbol  error 
causes  only  one  bit  error  as  follows: 


P  P 

s  _  s 


log?  M 


(0.31) 


where  k  is  the  number  of  bits  and  M  =  2k  is  the  size  of  the  constellation.  For  QPSK  the 
symbol  error  rate  is  given  by  [21] 
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For  higher  order  PSK  modulation  Ps  can  be  approximated  by 
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where  E  =  Eh  (log?  M)  is  the  energy  per  symbol.  As  seen  from  Equation  4.4,  the  bit  er¬ 
ror  rate  can  always  be  derived  by  dividing  Ps  by  the  number  of  bits. 


An  OFDM  system  using  QPSK  and  the  parameters  listed  in  Table  6  was  simu¬ 
lated  under  AWGN  conditions.  The  received  QPSK  constellation  is  shown  in  Figure  20 
for  a  a  value  of  0.02. 
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#  of  car¬ 
riers 

#  of  sym¬ 
bols 

Channel 

model 

seed 

Interlever 

matrix 

a  range 

48 

20000 

1 

33 

[144  139] 

0-0.08 

Table  6.  AWGN  Channel  Simulation  Parameters 
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Figure  20.  Effect  of  AWGN  on  Received  Signal  Constellation  for  a  =  0.02  (Before 

Differential  Decoding) 

Figure  21  shows  the  simulation  BER  plot  under  AWGN  conditions.  Clearly,  the 
errors  decrease  as  the  Eb/N0  increased.  This  result  is  used  in  later  simulations  as  a  refer¬ 
ence  for  comparing  the  results  with  frequency  errors  (see  Figures  29,  30,  34  and  35,  dis¬ 
cussed  later  in  the  thesis). 
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Figure  21 .  BER  versus  Eb/Na  Plot  for  a  QPSK  Simulation  in  AWGN  Channel 

3.  OFDM  Performance  in  Mobile  Channel  (Channel  Model  2) 
a.  Mobile  Channels 

Mobile  channels  are  composed  of  both  multipath  and  AWGN  components. 

The  multipath  component  of  mobile  Channels  1  and  2  has  a  delay  line  with  8  taps.  The 
input  time  domain  signal  is  filtered  by  4  parallel  7th  order  FIR  filters,  whose  coefficients 
are  fixed  according  to  [28].  These  filters  implemented  in  m-file  cvdd  as  in  Figure  22.  The 
FIR  filter  output  is  computed  using 

yi{n)  =  ^Jaikx{n-k)  i  =  1,2, 3, 4  (0.34) 

k= 0 

where  aik  are  the  filter  coefficients  and  x(n)  are  the  time- domain  input  samples.  The 
filter  coefficients  used  for  each  of  the  four  FIR  filters  are  fisted  in  Table  7. 
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k 

0 

1 

2 

3 

4 

5 

6 

7 

«i  k 

-0.013824 

0.054062 

-0.157959 

0.616394 

0.616394 

-0.157959 

0.054062 

0.013824 

a2k 

0.003143 

-0.019287 

0.1008 

-1.226364 

1.226364 

-0.1008 

0.019287 

-0.003143 

a3k 

0.055298 

-0.216248 

0.631836 

-0.465576 

-0.465576 

0.631836 

-0.216248 

0.055298 

a4  k 

-0.012573 

0.077148 

-0.403198 

0.905457 

-0.905457 

0.403198 

-0.077148 

0.012573 

Table  7.  FIR  Filter  Coefficients  for  Mobile  Channels  1  and  2 

*(«) 


Figure  22.  Modeling  of  Mobile  Channels  1  and  2  Using  FIR  Filters  and  Delay  Ele¬ 
ments  D. 

Multipath  components  of  mobile  channels-3  and  4  use  6-tap  FIR  filters  but 
they  have  time-varying  coefficients  with  Jakes  spectrum  [29].  The  filtered  output  is  de¬ 
rived  simply  by  summing  the  delayed  paths  with  the  direct  path. 

All  multipath  components  also  include  power  loss  of  received  signal  levels 
and  Doppler  frequency  shifting  characteristics.  (The  delay,  power  loss  and  Doppler  vec¬ 
tors  can  be  altered  to  reflect  the  desired  channel  characteristics.) 

b.  Mobile  Channel  Simulations: 

The  goal  of  this  simulation  step  was  to  see  the  mobile  channel  effects  on 
the  transmitted  symbols.  The  four  different  mobile  channel  models  were  simulated  and 
the  performance  curves  for  each  of  them  presented.  The  one  that  causes  the  most  degra¬ 
dation  was  chosen  to  be  able  to  simulate  the  characteristics  of  a  mobile  outdoor  environ¬ 


ment. 
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Mobile  channel  simulations  were  conducted  with  the  parameters  in  Table 
8.  Figure  23  shows  the  received  signal  constellations  after  each  mobile  channel  before 
differential  decoding  for  a  noise  variance  value  of  a  =  0.01 .  Since  the  mobile  channel 
causes  constructive  and  destructive  effects  on  the  OFDM  symbols,  the  received  constella¬ 
tion  points  are  scattered  from  their  transmitted  positions.  If  the  differential  encod¬ 
ing/decoding  were  not  implemented,  a  large  number  of  received  symbols  would  be  de¬ 
coded  in  error  since  many  points  are  scattered  into  adjacent  phase  sectors. 


#  of  carri¬ 
ers 

#  of  sym¬ 
bols 

Channel 

model 

seed 

Interlever 

matrix 

CT 

48 

20000 

21,22,23,24 

33 

[144  139] 

0.01 

Table  8.  Mobile  Channel,  Simulation  1  Parameters 
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(c)  (d) 

Figure  23.  Effect  of  Mobile  Channels  on  Received  Signal  Constellations  (Before  Dif¬ 
ferential  Decoding),  (a)  Mobile  Channel  1  (b)  Mobile  Channel  2  (c)  Mobile 
Channel  3  (d)  Mobile  Channel  4. 

Figure  24  shows  the  received  QPSK  constellations  after  differential  de¬ 
coding.  The  scattered  points  are  placed  within  their  previous  phase  sectors,  thereby  re¬ 
ducing  the  number  of  errors.  Figure  25  illustrates  the  effect  of  the  mobile  channels  on  the 
magnitude  plot  of  the  received  symbols.  As  previously  discussed,  the  sub-carriers  in  the 
middle  of  the  OFDM  spectrum  experience  more  multipath  loss  than  the  sub -carriers  on 
the  edges  of  the  spectmm  because  the  sub-carriers  in  the  middle  have  interfering  sub¬ 
carriers  on  both  sides.  The  effects  of  multipath  can  be  observed  in  the  form  of  magnitude 
distortion  in  Figure  25.  Also,  we  can  see  the  randomly  generated  peaks,  which  can  be  at¬ 
tributed  to  the  use  of  the  Jakes  spectmm  for  mobile  Channels  3  and  4. 
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(C)  (d) 

Figure  24.  Received  Signal  Constellation  Affected  by  Mobile  Channels  (After  Differ¬ 
ential  Decoding),  (a)  Mobile  Channel  1  (b)  Mobile  Channel  2  (c)  Mobile 
Channel  3  (d)  Mobile  Channel  4. 
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Figure  25 .  Effect  of  Mobile  Channels  on  the  Magnitude  Plot  of  Received  Symbols,  (a) 

Mobile  Channel  1  (b)  Mobile  Channel  2  (c)  Mobile  Channel  3  (d)  Mobile 

Channel  4. 

BER  plots  of  the  OFDM  system  under  the  mobile  channel  conditions  were 
obtained  by  configuring  the  system  with  the  parameters  in  Table  7.  A  range  of  a  values 
from  0  to  0.08  was  used  to  obtain  these  plots. 

Figure  26  shows  the  BER  versus  Eb /No  plot  of  each  mobile  Channel.  By 

observing  the  BER  curves,  it  can  be  seen  that  mobile  Channel  1  and  2  required  more  sig¬ 
nal  energy  relative  to  the  other  channel  models  for  the  same  BER.  Mobile  Channels  1  and 
2  had  approximately  the  same  performance,  and  we  randomly  chose  mobile  Channel  1  as 
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the  mobile  channel  model.  This  figure  justifies  that  Doppler  frequency  shift  have  negligi¬ 
ble  effects  on  the  performance  of  OFDM  systems  since  we  get  similar  BER  curves  for 
mobile  Channel  1  and  2.  The  Doppler  frequency  for  mobile  Channel  1  was  15  Hz  and  for 
mobile  Channel  2  it  was  5  Hz,  according  to  [29]. 


Figure  26.  BER  versus  Eb/Na  Plots  for  QPSK  Simulation  in  Mobile  Channels. 

4.  Performance  with  Frequency  offset 
a.  Noise  Free  Channel 

This  step  of  the  simulation  process  was  to  verify  if  it  were  possible  to  im¬ 
plement  the  frequency  offset.  As  previously  mentioned  in  Chapter  m,  for  the  frequency 
offset  case,  9  it)  is  deterministic  and  is  calculated  from27t  A F  t +0O ,  where  A F  is  the  car¬ 
rier  offset.  In  the  literature,  however,  frequency  offset  is  usually  expressed  as  a  ratio  of 
the  frequency  spacing  (AF/Rs) .  The  performance  degradation  caused  by  frequency  oflset 
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is  also  expressed  in  terms  of  this  relative  frequency  offset.  Therefore,  in  the  simulations, 
the  same  approach  was  used,  and  each  time  domain  sample  of  the  OFDM  symbols  was 
exposed  to  the  effective  value  of  the  relative  frequency  offset. 

To  see  the  frequency  offset  effects  exclusively,  a  noise- free  channel  was 
needed.  Thus,  Channel  0  was  chosen.  The  easiest  way  to  observe  the  effect  of  frequency 
offset  was  to  compare  the  signal  constellations  of  the  transmitted  and  received  symbols. 
The  code  to  simulate  a  channel  with  only  frequency  offset  was  configured  with  the  pa¬ 
rameters  in  Table  9.  A  relatively  small  frequency  offset  value  was  chosen  in  order  to  ob¬ 
serve  its  effect  clearly. 


#  of  car¬ 
riers 

#  of  sym¬ 
bols 

Channel 

model 

seed 

Interleaver 

matrix 

A  FIR, 
(%) 

48 

8000 

0 

22 

[167  48] 

0.1 

Table  9.  Frequency  Offset  in  Noise  Free  Channel  Simulation 


Hie  received  QPSK  signal  constellation  before  differential  decoding  is 
shown  in  Figure  27  with  the  frequency  offset.  As  seen  in  the  figure,  the  center  of  the  con¬ 
stellation  points  is  shifted  from  their  original  positions.  Figure  28  shows  that  the  differen¬ 
tial  decoding  algorithm  successfully  placed  the  scattered  points  within  their  previous  sec¬ 
tors.  Due  to  the  small  frequency  offset  value,  the  errors  were  easily  corrected.  However, 
for  wide  ranges  of  frequency  offsets,  as  the  frequency  offset  value  increased,  so  did  the 
errors. 
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Figure  27.  Received  Signal  Constellation  in  Channel  Model  0  with  Frequency  Offset  of 
0.1%  of  Frequency  Spacing  (Before  Differential  Decoding) 


Figure  28.  Received  Signal  Constellation  in  Channel  Model  0  with  Frequency  Offset  of 
0.1%  of  Frequency  Spacing  (After  Differential  Decoding) 
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b.  AWGN  Channel 

This  simulation  step  is  important  because  all  theoretical  expressions  for 
frequency  offset,  in  the  literature,  were  derived  under  AWGN  channel  conditions.  Thus, 
the  goal  was  to  study  how  the  system’s  performance  deviated  from  the  theoretical  per¬ 
formance  with  increasing  frequency  offset. 

The  OFDM  system  was  configured  with  the  parameters  given  in  Table  10. 


#  of  car¬ 
riers 

#  of  sym¬ 
bols 

Channel 

model 

seed 

Interleaver 

matrix 

o 

range 

A  F/Rs 
(%) 

48 

20000 

1 

33 

[144  139] 

0-0.8 

0-0.4 

Table  10.  Frequency  Offset  in  AWGN  Channel  Simulation 


Figure  29  show  the  BER  performance  for  relative  frequency  offset  values 
ranging  from  0  to  0.4%  in  0.1%  increments.  As  can  be  seen  from  the  figure,  above  0.4% 

frequency  offset  BER  reaches  1 0  1 ,  which  is  unacceptably  high. 


Figure  29.  BER  versus  Eh  /Na  Plots  for  QPSK  in  AWGN  Channel  with  Relative  Fre¬ 
quency  Offset  Values  from  0%  to  0.4% 
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In  Figure  29,  the  required  Eb/Na  values  increased  paraboloically  as  the 
relative  frequency  offset  values  were  increased;  we  compare  these  with  the  BER  plot  for 
AWGN  simulation.  This  was  an  expected  result  because,  by  looking  at  Equation  3.17,  the 
SNR  degradation  increased  with  the  square  of  the  relative  frequency  offset. 

c.  Mobile  Channel 

The  mobile  channel  models  were  intentionally  chosen  to  be  severe  chan¬ 
nels  to  represent  the  outdoor  mobile  environments.  In  this  step,  the  goal  was  to  see  the 
frequency  offset  effect  in  a  severe  mobile  channel.  Channel  1  is  selected  for  this  reason. 
The  system  was  configured  with  the  same  parameters  presented  in  Table  9  but  mobile 
Channel  1  was  used  for  this  simulation. 

Figure  30  shows  the  BER  plot  of  the  system  in  mobile  channel  1  for  rela¬ 
tive  frequency  offset  values  ranging  from  0  to  0.4%. 


Figure  30.  BER  versus  Eb/Na  Plots  for  QPSK  in  Mobile  Channel  1  with  Relative  Fre¬ 
quency  Oflset  Values  from  0%  to  0.4% 
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In  Figure  30  the  dotted  lines  are  from  a  frequency  offset  simulation  in  the 
AWGN  channel  and  the  solid  lines  are  from  a  frequency  offset  simulation  in  mobile 
channel  1.  We  can  observe  the  degradation  caused  by  the  mobile  channel  by  comparing 
the  dotted  and  the  solid  lines  and  also  the  OFDM  system’s  performance  in  mobile  Chan¬ 
nel  1  can  be  seen. 

5.  Performance  with  Phase  Noise 
a.  Noise-Free  Channel 

As  the  first  step,  the  noise- free  channel  was  used  to  study  the  phase  noise 
effect  exclusively.  As  mentioned  previously,  for  the  phase  noise  case,  0  (t)  was  assumed 

to  be  a  Wiener  process  with  ii{0  (r)]  =  0  and  £  j(0  (/;  +  /)  -0  [to ))  j  -  471  [1 1/| ,  where 

fj  (Hz)  is  the  one-sided  3-dB  line  width  of  the  Lorentzian  power  density  spectrum  of  the 
free- running  carrier  generator.  A  Wiener  process  was  created  in  Matlab  with  zero  mean 
and  variance  4n\i  |r|  with  |j  as  a  percentage  of  the  frequency  spacing.  The  first  simula¬ 
tion  was  run  with  the  configuration  parameters  given  in  Table  1 1.  Figure  31  shows  the 
generated  phase  noise  as  a  Wiener  Process. 


#of 

carriers 

#  of  sym¬ 
bols 

Channel 

model 

seed 

Interlever 

matrix 

P/R, 

<%) 

48 

8000 

0 

222 

[167  481 

0.01 

Table  1 1.  Phase  Noise  in  Noise  Free  Channel  Simulation 
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Figure  31.  Phase  Noise  Generated  as  a  Wiener  Process 

Two  different  phase-noise  sub-blocks  were  created.  The  first  began  the 
noise  generation  process  with  the  first  OFDM  symbol  and  ended  when  the  last  time  do¬ 
main  sample  of  the  last  OFDM  symbol  transmitted,  and  the  second  created  phase  noise 
for  each  OFDM  symbol  individually,  i.e.,  one  Wiener  process  for  each  80  time  domain 
samples.  The  first  method  was  used  during  the  simulations  since  it  more  accurately  re¬ 
flected  the  practical  operation  of  a  communication  system. 

Figure  32  shows  the  received  QPSK  signal  constellation  as  a  result  of 
phase  noise,  before  differential  decoding.  Figure  33  shows  the  corrected  signal  constella¬ 
tion  after  differential  decoding.  The  differential  decoding  algorithm  successfully  placed 
the  scattered  points  within  their  respective  sectors. 
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Figure  32.  Received  Signal  Constellation  in  Channel  Model  0  with  fj  /  R  =  0.0001 . 

(Before  Differential  Decoding) 


Figure  33.  Received  Signal  Constellation  in  Channel  Model  0  with  [)//?  =  0.0001 . 

(After  Differential  Decoding) 
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b.  AWGN  Channel 

To  study  the  phase  noise  effect  on  the  system,  the  code  was  configured 
with  the  parameters  in  Table  11.  Figure  34  shows  the  BER  plots  of  the  system  in  AWGN 
channel  for  normalized  bandwidth  ( (3  /  )  values  from  0.07%  to  0.28%  in  0.07%  incre¬ 

ments. 


#of 

carriers 

#of 

symbols 

Channel 

model 

seed 

Interleaver 

matrix 

a 

range 

p/«, 

<%> 

48 

20000 

1 

33 

fl44  1391 

0-0.8 

0-0.0028 

Table  12.  Phase  Noise  in  AWGN  Channel  Simulation 


Figure  34.  BER  versus  Eb/Na  Plots  for  QPSK  in  AWGN  Channel  with  Normalized 

Bandwidth  (3 / Rs  Values  from  0.0%  to  0.28%. 
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Figure  34  suggests  a  linear  increase  in  the  required  Eb  /N()  values  to  ob¬ 
tain  the  same  error  rates  as  the  normalized  bandwidth  (1  / Rs  was  increased.  This  was  an 
expected  result  because  from  Equation  3.14,  the  SNR  degradation  increases  linearly  with 
the  increasing  normalized  bandwidth. 

c.  Mobile  Channel 

This  step  was  performed  to  investigate  the  effects  of  phase  noise  in  a  se¬ 
vere  mobile  channel.  The  same  parameters  in  Table  12  were  used,  but  the  simulation  was 
performed  under  mobile  Channel  1. 

Figure  35  shows  the  BER  performance  of  the  system  in  mobile  Channel  1 
for  normalized  bandwidth  values  ranging  from  0  to  0.28%  in  0.07%  increments.  The  dot¬ 
ted  lines  are  for  the  phase  noise  simulation  in  the  AWGN  channel  and  the  solid  lines  are 
for  a  phase  noise  simulation  in  the  mobile  Channel  1.  From  Figure  35,  we  can  observe  the 
degradation  due  to  mobile  channel  by  comparing  the  phase  noise  simulation  results  in  an 
AWGN  channel  (dotted  lines)  with  the  phase  noise  simulation  results  in  the  mobile  chan¬ 
nel  (solid  lines).  The  OFDM  system  under  mobile  channel  conditions  requires  an  Eh/Na 
value  of  approximately  2  dB  more  than  the  AWGN  case. 
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Figure  35.  BER  versus  Eb/Na  Plots  for  QPSK  in  Mobile  Channel  1  with  Normalized 

Bandwidth  p / Rs  Values  from  0.0%  to  0.28%. 


C.  SUMMARY 

This  chapter  has  described  a  Matlab- based  OFDM  communication  system  that 
can  meet  the  bit  rate  and  performance  requirements  of  the  IEEE  802.1  la  standard.  After 
the  verification  of  the  code,  several  simulation  steps  were  performed  under  a  variety  of 
channel  conditions.  Some  of  these  simulations  especially  the  ones  in  steps  4  and  5  took 
seven  hours,  which  was  a  limiting  factor  in  the  simulation  process. 

The  simulation  results  did  not  exactly  match  the  theoretical  results;  however,  they 
reflected  the  expected  behavior  according  to  the  theory.  On  the  other  hand,  the  results 
were  similar  to  those  reported  in  studies  for  the  same  type  of  simulations. 
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V.  CONCLUSION 


The  objective  of  this  thesis  was  to  investigate  the  frequency  offset  and  phase 
noise  effects  on  OFDM  based  communication  systems.  This  objective  was  accomplished 
by  investigating  the  effects  of  frequency  offset  and  phase  noise  on  the  performance  of  the 
system  through  a  comprehensive  simulation  study. 

A.  SUMMARY  OF  THE  WORK  DONE 

A  comprehensive  background  of  OFDM  based  communication  systems  was  pre¬ 
sented,  including  the  sub-carrier  orthogonality  requirement  and  IFFT/FFT  processes.  A 
theoretical  analysis  of  frequency  errors  and  their  effects  on  overall  system  performance 
were  examined  based  on  the  work  reported  in  the  literature. 

A  Matlab  simulation  of  an  OFDM  communication  system  was  developed.  The 
simulation  studies  were  performed  for  the  following  channels:  ideal,  AWGN,  and  mobile. 
The  mobile  channel  includes  multipath  effects  in  AWGN.  The  effects  of  frequency  offset 
and  phase  noise  were  studied  for  all  channel  models  by  observing  the  deviation  from  the 
AWGN  case  with  increasing  frequency  offset  and  phase  noise.  Simulations  studied  the 
BER  performance  plots  of  OFDM  system  affected  by  frequency  errors  under  AWGN  and 
mobile  channel  conditions. 

B.  SIGNIFICANT  RESULTS  AND  CONCLUSIONS 

The  results  were  satisfactory  in  terms  of  the  OFDM  system  exhibiting  the  ex¬ 
pected  behavior.  However,  there  were  some  discrepancies  with  the  theoretical  results. 
These  discrepancies  stem  from  many  different  factors.  To  make  a  one-to-one  comparison 
of  results  from  different  studies,  the  simulations  must  be  performed  under  the  same  con¬ 
ditions.  This  was  not  possible  in  this  thesis  because  the  conditions  were  either  different  or 
not  known  in  detail.  Therefore,  observing  the  expected  trend,  rather  than  obtaining  the 
exact  numerical  results,  was  enough  to  realize  the  objectives  of  this  thesis. 

From  the  theoretical  analyses  given  in  Chapter  3,  we  have  seen  that  the  SNR  deg¬ 
radation  due  to  the  Doppler  shift  is  generally  negligible.  The  simulations  support  this 
conclusion.  In  Figure  26,  we  can  see  that  mobile  Channel  1  and  2  have  similar  BER 
curves  even  though  the  former  has  15  Hz  and  the  latter  has  5  Hz  of  Doppler  frequency. 
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The  theoretical  analyses  for  frequency  offset  showed  that  the  OFDM  system’s 
performance  degrades  with  the  square  of  the  relative  frequency  offset.  This  result  was 
observed  in  our  simulations.  From  Figure  29,  a  frequency  offset  of  0.1%  required  1  dB 
more  Eb/Na  value  than  the  AWGN  case,  and  a  frequency  offset  of  0.2%  required  2.8  dB 
more  Eb/Na  for  the  same  BER. 

In  the  phase  noise  case,  the  theoretical  analyses  showed  that  the  OFDM  perform¬ 
ance  degraded  linearly  as  we  increased  the  normalized  bandwidth.  The  simulation  result 
in  Figure  34  supports  this  point. 

Simulation  results  indicated  that  mobility  causes  performance  degradation.  Fig¬ 
ures  30  and  35  show  that  the  performance  of  OFDM  under  mobility  conditions  is  worse 
compared  to  that  of  an  AWGN  channel. 

C.  SUGGESTIONS  FOR  FUTURE  STUDIES 

This  thesis  was  conducted  under  a  certain  set  of  conditions.  The  effects  of  chang¬ 
ing  these  conditions  on  the  performance  of  OFDM  systems  open  a  wide  range  of  research 
topics. 

A  baseband  implementation  of  the  communication  system  was  simulated  in  this 
thesis.  A  future  study  may  focus  on  an  actual  RF  implementation  to  observe  the  system’s 
performance  in  a  more  realistic  scenario. 

Differential  encoding/decoding  technique  was  implemented  in  this  thesis.  This 
technique  may  be  replaced  with  coherent  detection  techniques  to  study  the  system’s  sen¬ 
sitivity  to  synchronization  issues  including  the  use  of  pilot  tones.  This  study  may  help 
evaluate  if  the  added  complexity  of  the  coherent  techniques  pays  off  in  terms  of  increased 
performance. 

This  work  focused  on  the  effects  of  frequency  errors  and  studied  only  four  multi- 
path  channel  models.  A  more  comprehensive  study  on  the  effects  of  mobile  channels  may 
be  considered  in  a  future  work.  Such  a  study  will  help  in  the  adoption  of  OFDM  for  the 
fourth- generation  mobile  communications. 
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QPSK  was  used  in  this  thesis.  The  effects  of  other  signaling  schemes  on  the  sys¬ 
tem’s  performance  may  be  studied  to  find  an  optimum  signaling  scheme.  Also  of  interest 
would  be  a  study  of  the  effects  of  changing  the  OFDM  system  parameters  used  in  this 
thesis,  especially  the  guard  interval  and  the  number  of  OFDM  sub -carriers  on  the  per¬ 
formance  of  the  system. 
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APPENDIX  A.  MATLAB  CODE  EXPLANATION 


This  appendix  explains  the  function  of  each  m-file  and  their  relations  with  each 
other.  The  code  is  separated  into  three  parts:  transmitter,  channel  and  receiver. 

COFDMSIM:  This  m-file  is  the  main  simulation- governing  file.  The  main  con¬ 
figuration  parameters,  such  as  sub -carrier  number,  FFT  point  number,  noise  variance 
range,  multi- path  losses,  multi-path  delays,  Doppler  frequencies,  guard  interval  length, 
channel  models,  seed  value,  interleaver  matrix  dimensions,  time/lfequency  differential 
encoding  and  constellation  sizes  are  chosen  here.  After  the  system  is  configured  with 
these  parameters,  the  m-file  CHANCDL  is  called. 

CHANCDL:  This  m-file  performs  the  simulations  depending  on  the  channel 
model  chosen  in  COFDMSIM.  This  m-file  calls  CDRCDLFT. 

A.  TRANSMITTER 

1.  CDRCDLFT 

This  is  the  COFDM  encoder  with  CDL  interleaver  and  is  called  to  compute  the 
frequency  array  of  prearranged  modulation  values  and  the  matrix  of  differentially  en¬ 
coded  complex  values.  It  calls  MLTPL. 

MLTPL  computes  the  possible  interleaver  matrix  dimensions.  CDRCDLFT  cre¬ 
ates  an  error  message  if  the  interleaver  matrix  cannot  accommodate  all  the  created  sym¬ 
bols. 

MARYMSG  generates  the  required  number  of  bits  randomly  and  then  sends  them 
to  the  convolutional  encoder.  MARYMSG  calls  CNV_ENCD  and  BM. 

CNV_ENCD  convolutionaly  encodes  the  randomly  generated  bits  in  MARYMSG. 
BM  then  converts  the  convolutionaly- encoded  bits  to  m-ary  symbols.  ( m  value  was  cho¬ 
sen  in  COFDM.) 

As  the  last  step,  MARYMSG  creates  the  m-ary  msg  matrix,  which  contains  the 
message  symbols.  After  obtaining  the  message  symbols  from  MARYMSG,  CDRCDLFT 
calls  CDLILV,  which  interleaves  the  message  symbols.  During  the  interleaving  opera¬ 
tion,  ROTM  is  called  to  rotate  the  input  vector. 
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With  the  creation  of  the  interleaved  message  matrix,  CDRCDLFT  calls  BM  again 
and  creates  n-ary  message  matrix  that  are  to  be  transmitted.  (The  choice  of  n  was  also 
made  in  COFDM  and  it  determines  the  constellation.  In  the  simulations,  n  =  2  for  QPSK 
signaling  was  used.) 

DIFCDRFT  performs  time  or  frequency  differential  encoding,  as  required.  If  time 
differential  encoding  is  chosen,  a  cumulative  summation  on  each  column  is  executed.  For 
frequency  differential  encoding,  a  cumulative  summation  on  each  row  is  performed.  The 
cumulative  summation  is  executed  such  that  the  first  element  is  added  to  the  next  element 
and  replaces  the  second  element  in  the  current  summation,  then  adds  the  current  sum  to 
the  third  element  and  replaces  it  also  with  the  current  sum.  All  the  sums  are  modulo -A. 

After  deriving  the  differentially  encoded  matrix,  its  elements  are  encoded  as  A- ary 
complex  modulation  values.  The  unit  circle  is  divided  into  the  number  phase  partitions 
depending  on  the  signaling  constellation  and  each  element  of  the  differentially  encoded 
matrix  is  mapped  to  one  of  these  phases.  A  reference  row  of  ones  is  appended  to  the  mes¬ 
sage  array  to  provide  a  reference  starting  point  for  the  receiver,  in  the  case  of  time  differ¬ 
ential  encoding,  and  two  reference  columns  of  ones  are  inserted  in  the  case  of  frequency 
differential  encoding.  These  reference  symbols  are  stripped  off  in  the  receiver  during  dif¬ 
ferential  decoding. 

CMV2FA  is  the  last  m-file  CDRCDLFT  calls. 

CMV2FA  rearranges  the  differentially  encoded  complex  modulation  values  into  a 
frequency  array  to  prepare  them  for  OFDM  frequency  generation  by  the  IFFT. 

This  final  step  of  the  encoder  block  has  the  convolutionaly  encoded,  block  inter¬ 
leaved,  differentially  encoded  and  prearranged  frequency  array.  The  frequency  array  is 
returned  to  the  CHANCDL  m-file  for  the  IFFT  process.  TDA  is  the  m-file  responsible  for 
the  IFFT  process. 
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2.  TDA 

Performs  the  IFFT  processing  and  generates  the  OFDM  frequencies.  The  trans¬ 
mitted  OFDM  symbols  for  multi-path  effects  are  prepared  by  appending  a  periodic  guard 
interval  with  length  Ng.  ( Ng  is  chosen  in  COFDM  and,  guard  interval  details  were  dis¬ 
cussed  in  Chapter  II.) 

Hie  time  domain  samples  representing  the  OFDM  symbol  appended  with  the 
guard  interval  are  returned  to  CHANCDL  as  the  signal  to  be  transmitted  through  the 
channel.  Figure  36  shows  the  relations  of  m- files  in  the  transmitter  block  diagram. 


Figure  36.  Hierarchical  m-file  Block  Diagram  for  the  Transmitter 


B.  CHANNEL 

One  of  the  three  main  channel  models  can  be  chosen  and  combined  with  fre¬ 
quency  offset  or  phase  noise.  The  phase  noise  can  be  applied  prior  to  the  channel  to  simu¬ 
late  the  phase  noise  case. 

1.  Phase  Noise  Simulations 

The  time  domain  OFDM  symbol  samples  from  the  TDA  m-file  are  taken  and  sent 
to  the  phase  noise  m- files  of  PHSNS  or  phsns2. 

phsns2  generates  phase  noise  as  a  Wiener  process  and  then  applies  it  to  all  the 
samples  that  are  being  transmitted.  Then,  CHANCDL  sends  the  time  domain  samples  to 
the  desired  channel  model.  Figure  37  illustrates  the  place  of  the  phase  noise  block  in  the 
m-file  hierarchy. 
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Figure  37.  Phase  Noise  Creating  Block  in  the  m-file  Hierarchy 

2.  Channel  Model  Simulations 

Time  domain  samples  are  next  sent  to  the  channel  model  of  choice. 
a.  Channel  Model  0 

This  is  the  ideal  noise- free  channel  used  for  code  check  simulation  pur¬ 
poses.  Here,  the  output  of  the  transmitter  block  is  sent  directly  to  the  receiver.  This  model 
can  also  be  used  to  check  whether  frequency  offset  or  phase  noise  can  be  implemented. 
These  two  possibilities  are  illustrated  in  Figure  38. 


Figure  38.  Channel  Model  0  m-file  Hierarchy 
b.  Channel  Model  1 

If  channel  model  1  is  chosen,  then  CHANCDL  directs  the  time  domain 
samples  to  AWGN. 

AWGN  generates  a  matrix  of  complex  random  numbers  chosen  from  a 
normal  distribution  with  mean  0  and  variance  1.  This  matrix  has  the  same  dimensions  as 
the  time  domain  samples  matrix.  This  matrix  is  then  multiplied  by  the  noise  variance  pa¬ 
rameter  chosen  in  COFDM  and  added  to  the  time  domain  signal  matrix.  This  channel 
model  can  either  be  used  individually  or  as  a  component  of  a  mobile  channel. 
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c.  Mobile  Channel  Models 

Four  different  mobile  channels  are  created  by  combining  an  AWGN  chan¬ 
nel  with  multi- path  effects.  If  desired,  the  multi- path  channels  can  also  be  used  individu¬ 
ally  for  simulation  purposes.  The  choice  of  a  mobile  channel  model  is  determined  in  the 
COFDM  m-ffle. 

Mobile  Channel  1  and  2:  If  one  of  these  mobile  channels  is  chosen  in 
COFDM,  then  CHANCDL  sends  the  time  domain  signal  to  CHUHF.  CHUHF  functions 
as  the  multi- path  component  of  mobile  Channels  1  and  2.  It  calls  PI  JNF  to  realize  a  de¬ 
lay  line  and  set  up  the  multiple  delayed  paths. 

CVDD  creates  a  continuous  digital  delay  element  in  which  the  time  do¬ 
main  samples  are  filtered  using  an  eight-tap  FIR  filter  whose  tap  coefficients  are  a  func¬ 
tion  of  the  desired  delay. 

After  receiving  the  delayed  paths  from  DLINE,  CHUHF  next  calls  OFST, 
which  applies  a  frequency  offset  as  a  ratio  of  the  max  Doppler  shift  to  the  direct  path  of 
time  domain  samples. 

RAY _DOP  calculates  the  max  Doppler  shift  frequency  as  a  ratio  of  fre¬ 
quency  spacing.  A  sequence  of  complex  numbers  with  zero  mean  and  variance  0.5  is  cre¬ 
ated  having  a  Rayleigh  envelope  with  a  mean  square  value  of  1.  The  real  and  imaginary 
parts  are  generated  independently. 

Finally,  CHUHF  includes  the  power  losses  for  the  individual  multi-paths 
by  multiplying  each  loss  amount  by  the  respective  delay  line  output  arrays.  At  this  point, 
CHANCDL  has  the  time  domain  representation  of  the  transmitted  signal  plus  the  multi- 
path  effects.  The  time  domain  samples  are  then  sent  to  AWGN  as  the  second  component 
of  the  mobile  channel. 

Mobile  Channel  3:  If  this  model  is  chosen  in  COFDM,  then  CHANCDL 
routes  the  signal  to  channel_a,  which  models  different  Doppler  frequencies,  power  losses, 
envelope  (Ricean)  and  delays.  It  calls  jakes,  which  implements  FIR  filtering  by  using  a 
Hamming  window. 
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Mobile  Channel  4\  If  this  model  is  chosen,  CHANCDL  sends  our  signal  to 
channel_b,  which  is  the  last  mobile  channel  model  with  very  similar  characteristics  to 
model  3.  It  also  calls  jokes  for  filtering  purposes. 

Mobile  Channels  3  and  4  also  calls  AWGN  as  the  second  component  of 
the  mobile  channel. 

Figure  39  shows  hierarchical  architecture  of  all  the  mobile  channels.  For 
each  simulation,  it  is  only  possible  to  choose  one  mobile  channel.  However,  as  the  white 
blocks  suggests,  it  is  possible  to  combine  mobile  channels  with  phase  noise  or  frequency 
offset.  Mobile  channels  can  also  be  used  individually. 


Figure  39.  Mobile  Channels  m-file  Hierarchy 
Figure  40  summarizes  the  entire  channel  model  m-file  architecture. 
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Figure  40.  Summary  of  Channel  Model  m-file  Architecture 
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3.  Frequency  Offset  Simulations 

If  frequency  offset  simulations  are  desired,  CHANCDL  directs  the  time  domain 
samples  to  frofst,  after  the  channel  model  of  choice. 


frofst  applies  the  frequency  offset  to  each  of  the  transmitted  time  domain  samples 
as  a  percentage  of  the  frequency  spacing.  Figure  41  shows  the  place  of  frequency  offset 
generating  block  in  the  m-file  hierarchy. 


Figure  41 .  Frequency  Offset  Block  in  the  m-file  Hierarchy 


Finally,  all  the  m- files  in  the  channel  part  have  been  discussed  .The  last  part  of  the 
code  implements  the  receiver. 

C.  RECEIVER 

The  first  m-file  in  the  receiver  part  is  the  ITDA. 

IT  DA  takes  the  time  domain  samples,  removes  the  guard  interval  and  performs 
FFT  to  create  the  frequency  domain  samples.  The  decoding  functions  are  performed 
within  the  DECDRCDL. 

DECDRCDL  performs  the  decoding  and  deinterleaving  functions.  It  calls 
FA2CMA,  which  reconstmcts  the  frequency  array  as  a  complex  modulation  array  with 
the  correct  ordering  of  frequencies. DECDRCDL  next  calls  DFDCDRFT  to  decode  the 
complex  modulation  values  differentially. 

DFDCDRFT  takes  the  complex  modulation  values  from  FA2CMA  and  performs 
inverse  mapping  of  the  complex  numbers  to  N-  ary  symbols  based  on  the  value  chosen  in 
COFDM.  Differential  decoding  is  the  inverse  of  encoding  performed  in  the  transmitter. 

After  receiving  A- ary  message  from  DFDCDRFT,  DECDRCDL  calls  MB  and 
BM  respectively  to  reformat  the  received  symbols  into  m- ary  symbols. 
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The  next  step  in  the  receiver  block  is  to  deinterleave  the  message  symbol  array. 
DECDRCDL  calls  CDLDLV,  which  creates  the  deinterleaved  message  matrix  by  using 
the  ROTM  as  described  in  the  transmitter  block. 

These  deinterleaved  message  symbols  are  next  changed  to  binary  and  sent  to  the 
Viterbi  decoder. 


VITERBI  performs  the  Viterbi  decoding  and  during  the  process  it  calls  BIN2DECI 
and  DECI2BIN  m- files. 


Finally,  DECDRCDL  takes  the  received  bit  sequence  from  VITERBI,  reformats  it 
as  an  m-ary  message  sequence  and  creates  the  received  message  matrix.  Figure  42  shows 
the  receiver  m-file  block  diagram. 


Figure  42.  Hierarchical  m-file  Block  Diagram  for  the  Receiver 


At  this  step,  CHANCDL  successfully  governed  the  transmitter  and  receiver  op¬ 
erations.  Further  steps  are  taken  to  plot  the  necessary  data  to  be  used  for  analysis  pur¬ 
poses.  The  simulations  end  after  providing  the  simulations’  performance  curves. 


Figure  43  summarizes  the  entire  m-file  algorithm. 


Figure  43 .  General  m-  file  Block  Diagram 
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APPENDIX  B.  MATLAB  CODE 


In  this  appendix,  the  m- files  used  in  the  simulations  are  presented.  Descriptions  of 
the  m- files  included  here  can  be  found  in  Appendix  A. 

%AWGN 

% - 

%Title  :  Additive  White  Gaussian  Noise  (Channel  Model  1) 

%Author  :  Prof.  Paul  H.  Moose,  Naval  Postgraduate  School 
%Modified  by  :  Ahmet  Yasin  ERDOGAN,  Naval  Postgraduate  School 

% 

%INPUTS  : 

%  X  :  Time  domain  samples  of  the  transmitted  signal 
%  s  :  Seed  parameter 

%  sigma  :  AWGN  Noise  variance  parameter  for  calculating  Eb/No 

% 

%OUTPUTS: 

%  Y  :  Time  domain  samples  of  the  transmitted  signal  plus  AWGN 

% 

function  Y  =  awgn(X,  sigma,  s) 

% 

%Find  dimensions  of  the  input  array 
[rr,cc]=size(X); 

% 

randn('seed',s+30); 

%Generate  a  random  real  part 
wreal=randn(rr,cc) ; 

%Generate  a  random  imaginary  part 
randn('seed',s+40); 
wimg=i*randn(rr,cc) ; 

% 

%An  array  of  random  complex  entries  chosen  from  a  normal  distribution  with 
%mean  0.0  and  variance  1.0.  Array  dimensions  are  the  same  as  X. 

W=wreal+wimg; 

% 

%Random  noise  multiplied  by  the  sigma  factor  and  added  to  the  signal. 

Y =X+(sigma.  *W) ; 
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%BIN2DECI 

%- . . . . . - . 

% 

%Title  :  Binary  To  Decimal  Conversion 

%Author  :  Tan  Kok  Chye,  Naval  Postgraduate  School 

% 

%INPUT  : 

%  v_x  :  Binary  input 

% 

%OUTPUT: 

%  v_y  :  Decimal  output 

% 

function  v_y=bin2dcci(v_x) 
v_l=length(v_x); 
v_y=(v_l-l:-l:0); 
v_y=2.Av_y; 
v  y=v  x*v  y'; 
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%BM 

% 

% 

%Title  :  Binary  to  M-ary  Converter 

%Author  :  Prof.  Paul  H.  Moose,  Naval  Postgraduate  School 

%Modified  by  :  Ahmet  Yasin  ERDOGAN,  Naval  Postgraduate  School 

% 

% 

%INPUTS: 

%  q  :  Base  2  exponent  for  M-ary  symbol  generation 
%  v  :  Binary  data  vector 

% 

%OUTPUTS: 

%  m  :  M-ary  output  vector  in  decimal  notation 

% 

function  m=bm(q,v) 

% 

%Find  the  length  of  input  vector, v, and  determine  if  there  is  a  remainder 

%after  dividing  by  q 

n=length(v); 

r=rcm(n,q); 

% 

%If  there  is  no  remainder, don't  pad  v  input  vector.  Otherwise  add  the  appropriate 
%number  of  zeros  to  generate  a  code  word  with  an  exact  multiple  of  q  bits. 

% 

if  r==0 

v=v; 

else 

v=[v  zeros(l,q-r)]; 
end 

% 

%Place  least  significant  bit  of  the  symbol  on  the  left  end. 
map=l; 
for  j=l:q-l 
map=[map  2Aj]; 
end 
% 

%Remove  q  bits  at  a  time  from  v  to  generate  m-ary  symbol  values. 

n=length(v); 

p=round(n/q); 

A=zeros(q,p); 


A(:)=v; 

m=map*A; 

m_aiy_msg=m; 
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%CDLDLV 

% 

% 

%Title  :  CDL  Block  Deinterleaver 

%Author  :  Prof.  Paul  H.  Moose,  Naval  Postgraduate  School 
%Modified  by  :  Ahmet  Yasin  ERDOGAN,  Naval  Postgraduate  School 

% 

%— - 

%INPUTS: 

%  1  :  Intermediate  matrix  row  number 

%  k  :  Intermediate  matrix  column  number 

%  dcase:  Indicates  the  deinterleaving  method  to  be 

%  used  (9  different  cases).  In  our  thesis  only  block  interleaving  was  used. 

%  si  :  Input  message  string  to  be  deinterleaved 

%  SYNC  :  Frame  synchronization  bits  (Not  used  in  COFDM  simulation) 

% 

%OUTPUTS: 

%  s  :  Interleaved  output  string 

% 

function  s=cdldlv(l,k, dcase, si, SYNC) 

si(length(si)+l-length(SYNC):length(si))=zeros(l,length(SYNC)); 

N=length(si); 

if  l*k==N 

x=zeros(l,k); 

x(:)=si; 

K=(l:k)-1; 

CR=K.*(K+l)/2; 

L=(l:l)-1; 

RR=L.*(L+l)/2; 

% 

if  dcase==l 
for  kk=l:k 

x(:,kk)=rotm(x(:,kk),CR(kk)); 

end 

elseif  dcase==2 
for  kk=l:k 

[z,x( :  ,kk)]  =rotm(x( :  ,kk)  ,CR(kk)) ; 
end 

elseif  dcase==3 
for  kk=l:l 

x(kk, :  )=rotm(x(kk, :  ),RR(kk)) ; 
end 

elseif  dcase==4 
for  kk=l:l 

[z,x(kk, : )]  =rotm(x(kk, : )  ,RR(kk) ) ; 
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end 

elseif  dcase==5 
for  kk=l:k 

x(:,kk)=rotm(x(:,kk),CR(kk)); 

end 

for  11=1:1 

x(ll,:)=rotm(x(ll,:),RR(ll)); 

end 

elseif  dcase==6 
for  kk=l:k 

[z,x(:,ld<:)]=rotm(x(:,kk),CR(kk)); 

end 

for  11=1:1 

x(ll,:)=rotm(x(ll,:),RR(ll)); 

end 

elseif  dcase==7 
for  kk=l:k 

x(:,ldc)=rotm(x(:,kk),CR(kk)); 

end 

for  11=1:1 

[z,x(ll,:)j=rotm(x(ll,:),RR(ll)); 

end 

elseif  dcase==8 
for  kk=l:k 

[z,x(:,ld<:)]=rotm(x(:,kk),CR(kk)); 

end 

for  11=1:1 

[z,x(ll,:)]=rotm(x(ll,:),RR(ll)); 

end 

end 

x=x'; 

s=x(:); 

s=s'; 

end 
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%CDLILV 

% 

% 

%Title  :  CDL  Block  Interleaver 

%Author  :  Prof.  Paul  H.  Moose,  Naval  Postgraduate  School 
%Modified  by  :  Ahmet  Yasin  ERDOGAN,  Naval  Postgraduate  School 

% 

%— - 

%INPUTS: 

%  1  :  Intermediate  matrix  row  number 

%  k  :  Intermediate  matrix  column  number 

%  dcase:  Indicates  the  deinterleaving  method  to  be 

%  used  (9  different  cases).  In  our  thesis  only  block  interleaving  was  used. 

%  si  :  Input  message  string  to  be  deinterleaved 

%  SYNC  :  Frame  synchronization  bits  (Not  used  in  COFDM  simulation) 

% 

%OUTPUTS: 

%  si :  Interleaved  output  string 

% 

%Subroutines  Used :  rotm.m 

%  - 

function  si  =  cdlilv(l,k, dcase, s, SYNC) 

N=length(s); 

if  l*k==N 

x=zeros(l,k); 

x=x'; 

x(:)=s; 

x=x'; 

Intermediate_mx=x; 

K=(l:k)-1; 

CR=K.*(K+l)/2; 

L=(l:l)-1; 

RR=L.*(L+l)/2; 

% 

if  dcase==l 
for  kk=l:k 

[z,x(:,kk)]=rotm(x(:,kk),CR(kk)); 

end 

elseif  dcase==2 
for  kk=  1  :k 

x(:,kk)=rotm(x(:,kk),CR(kk)); 

end 

elseif  dcase==3 
for  kk=l:l 

[z,x(kk, : )]  =rotm(x(kk, : )  ,RR(kk) ) ; 
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end 

elseif  dcase==4 
for  kk=l:l 

x(kk, :  )=rotm(x(kk, :  ),RR(kk)) ; 
end 

elseif  dcase==5 
for  11=1:1 

[z,x(U,:)]=rotm(x(ll,:),RR(ll)); 

end 

for  kk=l:k 

[z,x(:,kk)]=rotm(x(:,kk),CR(kk)); 

end 

elseif  dcase==6 
for  11=1:1 

[z,x(ll,:)]=rotm(x(ll,:),RR(ll)); 

end 

for  kk=l:k 

x(:,ldc)=rotm(x(:,kk),CR(kk)); 

end 

elseif  dcase==7 
for  11=1:1 

x(ll,:)=rotm(x(ll,:),RR(ll)); 

end 

for  kk=l:k 

[z,x(:,ld<:)]=rotm(x(:,kk),CR(kk)); 

end 

elseif  dcase==8 
for  11=1:1 

x(ll,:)=rotm(x(ll,:),RR(ll)); 

end 

for  kk=l:k 

x(:,kk)=rotm(x(:,kk),CR(kk)); 

end 

end 

x; 

si=x(:); 

si=si'; 

end 

si(length(si)-length(SYNC)+l:length(si))=SYNC; 
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%CDRCDLFT 

% 

% 

%Title  :  COFDM  Encoder  with  CDL  Interleaver 
%Author  :  Prof.  Paul  H.  Moose,  Naval  Postgraduate  School 
%Modified  by  :  Ahmet  Yasin  ERDOGAN,  Naval  Postgraduate  School 

% 

% 

%INPUTS: 

%  pic  :  Parameter  indicating  the  figure  number 
%  s  :  Seed  parameter 

%  freqno  :  Number  of  OFDM  frequencies  (sub -carriers)  used  in  each  message  array 
%  rintlv  :  Parameter  for  intermediate  matrix  row  number 
%  cintlv  :  Parameter  for  intermediate  matrix  column  number 
%  N  :  Number  of  FFT  frequency  sample  points, must  be  larger  than  freqno 
%  mary  :  Initial  M-ary  symbol  format  (OFDM  symbol  bit  length) 

%  nary  :  Final  N-ary  symbol  format  (PSK  symbol  bit  length) 

%  fort  :  Selects  either  frequency  (fort=l)  or  time  (fort=0)  differential  encoding 
% 

%OUTPUTS: 

%  Fa  :  Frequency  array  of  prearranged  modulation  values 
%  MD  :  Matrix  of  differentially  encoded  complex  values 
%(unity  magnitude  and  one  of  N-ary  possible  phases  (N-PSK)) 

%  B  :  Matrix  of  8- ary  symbols 

%  nsymno:  Number  of  N-ary  generated  symbols 

% 

%Subroutines  Used :  marymsgm,cdlilv.m,mb.m,bmm,difcdrft.m,cmv2fa.m 

% 

func- 

tion[Fa,MD,B_ce,B_random,nsymno]=cdrcdlft(pic,dcase,s, freqno, rintlv  ,cintlv,N,mary,na 
ry,fort); 

% 

%  Determine  whether  the  number  of  OFDM  frequencies  are  even,  indicated 
%  by  the  "freqno"  parameter.  If  odd,  go  to  error  message.  Odd  frequencies  are  not  al¬ 
lowed 

%  since  the  formation  of  the  frequency  array  is  symmetrical  and  even. 

% 

if  rem(freqno,2)~=0 

disp(FRROR:  The  number  of  matrix  columns,  representing  OFDM  frequencies,  must  be 
an  even  number!') 
elseif  rem(freqno,2)==0 

% 

%  Check  if  interleaver  matrix  dimensions  are  greater  than  freqno. 

%  If  not,  then  display  error  message  and  stop. 

% 
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if  (rintlv*cintlv)<(freqno) 
disp(") 

disp('ERROR:  The  row  and  column  interleave  parameters  are  not  compatible  with  #  of 
OFDM  frequencies!') 
disp(") 
else 

%  Calculate  the  row  symbol  number 
symno=rintlv*cintlv/freqno; 

% 

%  Display  error  message  if  symno  and  freqno  not  compatible  with  rintlv  and  cintlv  and 
stop. 

%  If  not  compatible, the  interleaver  function  does  not  work  correctly. 

% 

if  rem(symno ,  1 )  ~=0 
disp(") 

disp('ERROR:  The  row  and  column  interleave  parameters  are  not  compatible  with  # 
of  OFDM  frequencies!') 

disp('  For  the  entered  rintlv,  cintlv,  and  freqno  parameters,  the  calculated  symno  is:') 
disp(symno) 

multiesall=mltpl(rintlv, cintlv); 
multies=multiesall(  1  ,(2:length(multiesall)  - 1 )) ; 
disp('  Possible  choices  for  freqno  based  upon  rintlv  and  cintlv  are:') 
disp(") 
disp(multies) 
elseif  rem(symno,l)==0 
if  freqno  >=  N; 
disp(") 

disp('ERROR:  The  number  of  frequency  points,  N,  needs  to  be  increased  !') 

disp('N  must  be  larger  than:') 

disp(") 

disp(freqno) 

disp(") 

elseif  freqno  <  N; 

% 

%  Generate  a  random  message  matrix  of  m-ary  symbols,based  upon  parameter, mary. 
Nmbr_of_symbols  =symno*freqno; 

[B_ce,B_random]=marymsg(mary, symno, freqno); 
Rndm_m_ary_msg=B_random; 

% 

%  Perform  a  block  interleaving  function  on  the  matrix,  B,  with  rintlv  rows 
%  and  cintlv  columns. 

SYNC=[]; 

[Br  Bc]=size(B_ce); 

Bt=B_ce'; 

Bvect=Bt(:)'; 
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si=cdlilv(rintlv,cintlv,dcase,Bvect,SYNC); 

Bi=reshape(si,Bc,Br)' ; 

Lrtrlvd_array=Bi; 
ml  =bm(nary  ,mb(mary  ,Bi)) ; 
lengthml=length(ml); 
nsymno=lengthml ; 
remml=rem(lengthml  ,freqno); 
if  remml=0; 
ml=ml; 
else 

zero=zeros(freqno-remml); 

ml=[ml  zero(l,:)]; 

end 

Iength2ml=length(ml); 

m=(reshape(ml,freqno,length2ml/freqno))'; 

N_ary_msg=m; 

% 

%  Generate  a  differentially  encoded  matrix  of  complex 
%  values  with  unity  magnitude  and  one  of  (2An)  equal  phases. 
MDD=difcdrft(nary  ,m,fort) ; 

[MDm  MDn]=size(MDD); 

MD=MDD; 

Cmplx_mod_array=MDD; 

% 

%  Form  the  frequency  array  of  modulation  values  that  include  guard  interval. 
Fa=cmv2fa(N,MD) ; 

Freq_array=Fa; 

end 

end 

end 

end 
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%CHANCDL 

% 

% 

%Title  :  Simulations  for  AWGN  &  Multipath  Fading  Channel 
%Author  :  Dave  Roderick,  Naval  Postgraduate  School 
%Modified  By  :  Ahmet  Yasin  ERDOGAN,  Naval  Postgraduate  School 

% 

%Subroutines  Used :  cdrcdlft.m,tda.m,awgn.m,chuhf.m,itda.m,decdrcdl.m,check.m 

% 

function 

chancdl(chnmdl, wait, pic, dcase,s,freqno,rindv,cintlv,N,mary, nary  ,n,k,blklgth,Ng,sigs, loss 

,dly,dop,freqspace,fort) 

sigvect=sigs; 

klgth=length(k); 

chklp=l; 

symno=rintlv*cintlv/freqno; 

%%%%FOR  PHASE  NOISE%%%%% 

%for  perctB=0:  .0008:  .0032 

%%%%FOR  FREQUENCY  OFFSET%%%% 

%for  donci -0:0.0007 :  .0028 

errvect=[]; 

bervect=[]; 

freqerrmx=[]; 

errsperpr=[]; 

Es_No=[]; 

Eb_No=[] ; 
sermx=[]; 
bermx=[]; 
rowerrmx=[]; 

for  lp=l:length(sigvect); 

[xmt,modvals,B_ce,B_random,nsymno]=cdrcdlft(pic,dcase,s,freqno,rintlv,cintlv,N,mary, 
nary, fort); 

xmtifft=tda(Ng,xmt) ; 

%%%%%%%%%%%%FOR  PHASE  NOISE%%%% 
%sandnaphns=phsns2(xmtifft,freqspace,perctB); 

%xmtifft=sandnaphns ; 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

xmtptsM:length(xmtifft); 
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if  chnmdl=0 
sandn=xmtifft; 

elseif  chmndl=l 

disp(['Sigma=',num2str(sigvect(lp))]); 

sandn=awgn(xmtifft,sigvect(lp),s); 

elseif  chnmdl=21 

sandrnltpth=chuM(s+l,xmtifft,loss,dly,dop,N,freqspace); 

disp(['Sigma=',num2str(sigvect(lp))]); 

sandn=awgn(sandmltpth,sigvect(lp),s); 

elseif  chnmdl=22 

sandmltpth=chuhf(s+l  ,xmtifft,loss,dly,dop.N,freqspacc); 

disp(['Sigma=',num2str(sigvect(lp))]); 

sandn=awgn(sandmltpth,sigvect(lp),s); 

elseif  chmndl=23 

disp(['Sigma=',num2str(sigvect(lp))]); 
sandmltpth=channel_a(xmtifft) ; 
sandn=awgn(sandmltpth,sigvect(lp),s); 

elseif  chmndl=24 

disp(['Sigma=',num2str(sigvect(lp))]); 
sandm  1  tpth=chan  ncl_b(  x  ml  ifft) ; 
sandn=awgn(sandmltpth,sigvect(lp),s); 
end 

%%%%%%%%%FOR  FREQUENCY  OFFSET%%%%%% 

%sandnofst=frofst(sandn, doner) ; 

%sandn=sandnofst; 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

sandnfft=itda(Ng.  sandn) ; 

K=(length(modvals(  1 , :  )))/2; 

[rcvd,rcvd_hit,random_msg,random_hit.M.MM]=dccdrcdl(pic,dcasc,K,sandnfft,nsymno, 
freqno,rintlv,cintiv,rnary, nary, fort  ,B_random); 

%%Transmitted_msg=B_random; 

Transmitted_msg=random_msg ; 

Received_msg=rcvd; 

[er¬ 
rors, bit_error,freqerrs,errmx,rowerrs]=check(pic,random_msg,random_bit,rcvd,rcvd_bit, 
n,k(chklp),blklgth); 
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errvect=[errvect, errors] ; 
bervect=[bervect,bit_error] ; 
freqerrmx=[freqerrmx;freqerrs] ; 
rowerrmx=[rowerrmx;rowerrs] ; 

cmtEs_No=l/(2*N*(sigvect(lp)A2)); 
cmtEb_N  o=cmtEs_N  o ; 

%  %based  on  M=4  i.e.  for  coded  QPSK,  Eb=Es. 

Es_No=[Es_No,cmtEs_No] ; 

Eb_No=[Eb_No,cmtEb_No] ; 

Es_Nodb=  1 0*log  1 0(Es_N  o) ; 

Eb_Nodb=10*log  10(Eb_No) ; 
end 

sei-en'vcct/(symno*frcqno) ; 
ber=bervect/(2*symno*freqno) ; 
sermx=[sermx;ser] ; 
bermx=[bermx;ber] ; 
errsum=sum(errvect) ; 
errsperpr=[errsperpr,errsum] ; 
errmax=max(rowerrmx') ; 

% 

%  PLOTS 

% 

figure(pic+l) 
plofrmodvals,'*') 
hold  on; 
plot(0,0,'+') 
hold  off; 

%title( [Transmitted  Signal  Constellation', int2str(2Anary),  '-ary  (QPSK) ']) 

xlabel('Re'); 

ylabel('Im'); 

axis('square');orient  tall;grid 
pause(wait); 

%%%%Plot  of  FFT  POINTS%%%% 

%figure(pic+2) 

%plot([0:N  -  l],abs(xmt),'*') 

%title(['Frequency  Array  Plot  (number  of  FFT  frequency  points  are  ',int2str(N),')']) 
%xlabel(['Guard  interval  length  is  ',int2str(N-freqno)]) 

%axis('square');orient  tall;grid 
%pause(wait); 

% 

figure(pic+2) 
surf(abs(modvals)); 
shading  interp;grid;orient  tall 
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%title(  ['Magnitude  Plot  of  Transmitted  Signal']) 
xlabel('OFDM  Subcarrier  Number') 
ylabel('Symbol  Row  Number') 
zlabel('Signal  Magnitude') 
pause(wait); 

% 

figure(pic+3) 
plot(M,'*');hold  on; 
plot(0,0,'+');hold  off; 

%title(['Received  Signal  Constellation',  int2str(2Anary),'-ary  (QPSK),  (Before  Differential 

Decoding)']) 

xlabel('Re'); 

ylabel('Im'); 

orient  tall;axis('square');grid 
pause(wait); 

% 

figure(pic44) 
plot(MM,'+') 
hold  on; 
plot(0,0,'+') 
hold  off; 

%title(['Received  Signal  Constellation', int2str(2Anary), '-ary  (QPSK),  (After  Differential 

Decoding)']) 

xlabel('Re'); 

ylabel('Im'); 

orient  tall;axis('square');grid 
pause(wait); 

% 

figure(pic+5) 
surf(abs(M)); 
shading  interp 
grid;orient  tall 

%title(['Magnitude  Variation  of  Received  Signal  (Sigma=',num2str(sigvect(lp)),')']) 

xlabel('OFDM  Subcanier  Number') 

ylabel( 'Symbol  Row  Number') 

zlabel('Signal  Magnitude') 

pause(wait); 

% 

if  errsum~=0 

%%%%PERFORMANCE  PLOTS%%%% 

%theoretical  on  top  of  simulation  curves 
%q=length(Eb_No) 

%for  i=l:q 

%P=0.5*(erfc(Eb_No(l,i))*(l-0.25*erfc(Eb_No(l,i)))); 
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%Pb(l,i)=P; 

%end 

figure(pic+6) 

%semilogy(Eb_Nodb,Pb,'m');grid;hold  on; 
semilogy(Eb_Nodb,ber);hold  on;grid 
if  fort==l 

%tide('Performance  curve:  BER  vs.  Eb/No  for  QPSK  Signaling  (Freq.  Diff.Enc.)') 
elseif  fort==0 

%tide('Performance  curve:  BER  vs.  Eb/No  for  QPSK  Signaling  (Time.  DiffEnc.)') 
end 

xlabel('Eb/No(dB)'); 

ylabel('BER'); 

%legend('theorical  curvcfsimulation  curve') ;orient  tall 

% 

figure(pic+7) 

%semilogy(Es_Nodb,2*Pb,'m');hold  on 
semilogy(Es_Nodb,ser);grid;hold  on 
if  fort==l 

%title('Performance  curve:  SER  vs.  Es/No  (Freq.  Diff.Enc.)') 
elseif  fort==0 

%title('Performance  curve:  SER  vs.  Es/No  (Time.  Diff.Enc.)') 
end 

xlabel('Es/No(dB)'); 

ylabel('SER'); 

%legend('theoretical  curve','simulation  curve')  ;orient  tall 
end 
end 
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%channel_a 

% 

% 

%Title  :  Mobile  Channel  2 

%Author  :  Prof.Roberto  Cristi,  Naval  Postgraduate  School 
%Modified  By  :  Ahmet  Yasin  ERDOGAN,  Naval  Postgraduate  School 

% 

%Subroutines  Used :  jakes.m 

% 

function  yr=channcl_a(xt) 

[asd  asn]=size(xt); 

xt=xt.'; 

xt=xt(:).'; 

% 

%  yr=channel(xt) 

%  with  xt=transmitted  vector 
%  yr=received  vector 

%  length(yr)=length(xt) 

% 

%  if  the  program  gives  you  an  error,  just  increase  the  size  of  the 
%  transmitted  vector  "xt" 

Fs=l  .25*  10A6;  %  sampling  frequency  (Hz) 

Fd=200;  %  doppler  frequency  (Hz) 

%Parameters  for  C 
Td=[0.0,  0.25,0.5, 

PdB=[0.0,-1.0,-9. 

K=[0, 0,0, 0,0,0]; 

Td=Td*(10A(-6)); 
nd=round(Td*Fs); 

P=10.A(PdB/10); 

Np=length(xt); 

yr=zcros(sizc(xt)); 

for  k=  1  :length(Td) 
g=jakes(Fd,  Fs,  Np); 

s=sqrt(P(k)/(K(k)+ 1 )) ;  %  random  path 

m=sqrt(P(k)*K(k)/(K(k)+l));  %  direct  path 
total=s*g+m*exp(j  *2*pi*Fd/Fs*(0:length(xt)- 1 )) ; 
yr=yr+[zcros(  1  ,nd(k)) ,  xt(  1  :length(xt)- nd(k))] . *total; 
end 

yr=reshape(yr,asn,asd).'; 


%  time  delays  (in  microsec) 
%  Powers  (in  dB) 

%  Ricean  Factor 

%  time  delays  (in  sec) 

%  time  delay  in  samples 

%  Powers  (Finear) 


Channel  A 
,1.0, 1.9, 2.2]; 

0,- 10.0,- 15.0,-20.0]; 
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%channel_b 

% 

% 

%Title  :  Severe  Mobile  channel 

%Author  :  Prof.Roberto  Cristi,  Naval  Postgraduate  School 
%Modified  By  :  Ahmet  Yasin  ERDOGAN,  Naval  Postgraduate  School 

% 

%Subroutines  Used :  jakes.m 

% 

function  yr=channcl_b(xt) 

[asd  asn]=size(xt); 

xt=xt.'; 

xt=xt(:).'; 

% 

%  yr=channel(xt) 

%  with  xt=transmitted  vector 
%  yr=received  vector 

%  length(yr)=length(xt) 

% 

%  if  the  program  gives  you  an  error,  just  increase  the  size  of  the 
%  transmitted  vector  "xt" 

Fs=l  .25*  10A6;  %  sampling  frequency  (Hz) 

Fd=200;  %  doppler  frequency  (Hz) 

%  Parameters  for  Channel  B 

Td=[0.0,  0.25,9.0,13.0,  17.0,  20.0];  %  time  delays  (in  microsec) 

PdB=[-2.5,0.0,- 12.8,- 10,  -25.2,  -16];  %  Powers  (in  dB) 
K=[0.5,0.5,0,0,0,0];  %  Ricean  Factors 

Td=Td*(10A(-6));  %  time  delays  (in  sec) 

nd=round(Td*Fs);  %  time  delay  in  samples 

P=10.A(PdB/10);  %  Powers  (Finear) 

Np=length(xt); 

yr=zcros(size(xt)) ; 

for  k=  1  :length(Td) 
g=jakes(Fd,  Fs,  Np); 

s=sqrt(P(k)/(K(k)+ 1 )) ;  %  random  path 

m=sqrt(P(k)*K(k)/(K(k)-i-l));  %  direct  path 
total=s*g+m*exp(j  *2*pi*Fd/Fs*(0:length(xt)- 1 )) ; 
yr=yr+[  zeros!  1  ,nd(k)) ,  xt(  1  :length(xt)-  nd(k))] .  *total; 
end 

yr=reshape(yr,asn,asd).'; 
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%CHECK 

% 

% 

%Title  :  Source  and  Sink  Message  Checker 

%Author  :  Prof.  Paul  H.  Moose,  Naval  Postgraduate  School 

%Modified  by  :  Ahmet  Yasin  ERDOGAN,  Naval  Postgraduate  School 

% 

%-— 

func- 

tion[error_no,bit_error_total,freqerrs,errmx,rowerrs]=check(pic,x,xbit,y,ybit,n,k,blklgth) 

if  blklgth>n 

disp(") 

disp(ERROR!  The  block  length(blklgth),  must  be  equal  to  or  less  than  the  code  word 
length(n)') 

disp('Please  enter  a  smaller  value  for  blklgth,  or  change  n') 
disp(") 

elseif  blklgth<=n 
if  n<k 
disp(") 

disp('Error!  The  code  word  length(n)  must  be  equal  to  or  larger  than  the  information 
length(k)') 

disp('Please  enter  a  larger  value  for  n,  or  change  k  to  a  smaller  number.') 

disp(") 

elseif  n>=k 

First_matrix=x; 

Second_matrix=y; 

[rx  cx]=size(x); 

% 

%Compare  inputs  x  and  y  and  generate  error  matrix,  "errors" 
errors=(x~=y); 

First=xbit; 

Second=ybit; 

[rxl  cxl]=size(xbit); 

% 

%Compare  inputs  xbit  and  ybit  and  generate  BIT  error  matrix,  "bit_errors" 
bit_errors=(xbit~=ybit) ; 

% 

%Find  the  error  distribution  vs.  OFDM  frequencies 
freqerrs=sum(errors); 

% 

%Find  the  error  location  in  "errors"  where  element  in  x  and  y  differ. 
Error_locations=(find(errors))'; 

Error_number=sum(sum(errors)) ; 

Correct_symbl_num=(size(y,  1  )*size(y,2))-  Error_number; 

% 
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%Find  the  bit  error  location  in  "errors"  where  element  in  xl  and  yl  differ. 
%bit_Error_locations=(find(bit_errors))'; 

bit_error_total=sum(sum(bit_errors)); 

%Correct_bit_num=(size(y  1 , 1  )*size(y  1 ,2))-bit_Error_number; 

%Reed- Solomon  8- bit  symbol  correction  for  (n-k)/2  symbols 
symcon -floor)  (n-  k)/2) ; 
if  blklgth<=(n-k) 

disp('Error! !  !The  block  length  is  too  short  for  the  given  n  and  k  values') 
disp(") 

elseif  blklgth>(n-k) 
errtrans=errors'; 

% 

%Reshape  the  error  matrix  as  a  vector  of  errors 
errvect=errtrans( : ; 
blkrem=rem(length(errvect),blklgth); 
if  blkrem~=0; 

zeropad=zeros(blklgth-  blkrem) ; 
errvectpad=[errvect  zeropad(l,:)]; 
elseif  blkrem==0; 
errvectpad=errvect; 
end 
% 

blknos=length(errvectpad)/blklgth; 

errcorct=[]; 

errblksum=[]; 

% 

for  lp=l:blknos; 

errblk=errvectpad(((blklgth*(lp  - 1 ))+ 1 ) :  (blklgth*lp)) ; 
errblklgth=length(errblk); 

if  sum(errblk)<=symcorr; 
nocn-zcros(ciTblklgth) ; 
errblk=noerr(l,:); 

elseif  sum(en’blk)>symcorr; 
errblk=errblk; 
end 

errcorct=[errcorct  errblk]; 
errblksum=[errblksum  sum(errblk)]; 
end 

newerrvect=errcorct(  1  :length(errvect)); 
errtot=sum(newerrvect) ; 

RSerrs=(reshape(newerrvect,size(errors,2),size(errors,l)))'; 

% 

%Find  the  error  distribution  vs.  OFDM  Frequencies 
freqerrs=sum(RSerrs) ; 
errindex=(find(RSerrs))'; 
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RSerrtot=sum(errblksum); 

RSerrdif=Error_number-RSerrtot; 
errperblk=[(  1  :blknos)  ;errblksum] ; 

% 

%Check  if  x  and  y  are  the  same.  If  not,  display  error  message 
if  x=y; 

disp('There  are  no  errors! ! !') 

error_no=0; 

errmx=errors; 

rowerrs=sum(errors'); 

else 

disp('WARNING!:Errors  were  detected!1) 

disp(") 

if  n=k 

disp('WARNING!:  Since  n=k, there  is  no  R-S  error  correcting  possible') 

disp(") 

end 

disp(['For  the  given  input  parametersm- ,int2str(n),'and  k=',int2str(k),',the  Reed- 
Solomon  code  is  capable']) 

disp(['of  correcting  ',int2str(symcorr), 'errors.']) 
disp(") 

% 

%Check  if  xbit  and  ybit  are  the  same.  If  not,  display  error  message 
if  xbit==ybit; 

disp('There  are  no  bit  errors! ! !') 
bit_error_no=0; 
bit_errmx=bit_errors ; 
bit_rowerrs=sum(bit_errors'); 
else 

disp('WARNING!:Errors  were  detected!') 
disp(") 
if  n=k 

disp('WARNING!:  Since  n=k, there  is  no  R-S  error  correcting  possible') 
disp(") 
end 
end 

% 

%RS  code  was  able  to  correct  all  errors 
if  errtot==0 

Pre_RS_error_matrix=errors ; 

disp('EXCELLENT:  The  Reed- Solomon  code  corrected  all  detected  errors!') 

disp(['Originally  the  error  total  was:',int2str(Error_number)]) 

disp(") 

error_no=0; 

errmx=zeros(rx,cx) ; 
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rowerrs=sum(errmx') ; 


% 

%RS  code  was  able  to  correct  some  errors  but  not  all  of  them 
elseif  errtot<Error_number 
Pre_RS_error_matrix=error  s ; 

Post_RS_error_matrix=RSerrs ; 
errmx=RSerrs; 
rowerrs=sum(errmx') ; 

disp('The  Reed-Solemon  code  corrected  some  detected  errors,  but  not  all.') 

disp(['Originally  the  error  total  was  :  ',int2str(Error_number)]) 

disp(") 

disp([' After  R-S  decoding  ,  the  error  number  was  reduced  to:',int2str(RSerrtot)]) 
disp(") 

error_no=RSerrtot; 

disp(['Thc  total  number  of  correct  symbols  are:',int2str((size(y,l)*size(y,2))-RSerrtot)]) 
disp(") 

disp('The  error  number  distribution  per  block  number  is  :') 
disp(errperblk) 

% 

%RS  code  did  not  correct  any  errors 
elseif  errtot==Error_number 
Error_matrix=errors ; 
errmx=errors; 
rowerrs=sum(errors'); 

disp('The  Reed- Solomon  code  did  not  correct  any  errors.') 
disp('Perhaps  a  more  powerful  R-S  code  is  required.') 
disp(") 

disp(['The  total  number  of  error  occurrences  is:',int2str(Error_number)]) 
disp(") 

error_no=errtot; 

disp('The  error  number  distribution  per  block  number  is  :') 
disp(errperblk) 
end 
end 
end 
end 
end 

disp('. _ _'); 
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%CHUHF 

% 

%Title  :  UHF  Channel  Model  (multipath  Channel  Model  1) 
%Author  :  Prof.  Paul  H.  Moose,  Naval  Postgraduate  School 
%Modified  by  :  Ahmet  Yasin  ERDOGAN,  Naval  Postgraduate  School 

% 

function  y=chuhf(s,x,loss,dly,dop,N,freqspace) 

c=10.A(-loss./20); 

deltat=  1/(N  *freqspace) ; 

d=(dly .  *  .00000 1 )  ./deltat; 

e=dop  ./freqspace ; 

[L,Nt]=size(x); 

D=length(d); 

x=x.'; 

x=x(:).'; 

% 

%D  path  with  delays  from  d. 
xd=dline(x,d); 

[rr,cc]=size(xd); 

x=xd(l,:); 

% 

%  Offsets  direct  path  by  .7  of  max  doppler  freq. 
xo=of st(  .7  *e(  1 )  ,N,x) ; 

% 

%  First  path  with  no  fading, 
for  1=1  :D 

a=ray_dop(s  ,cc,N,e(l)); 

xd(l,:)=a.*xd(l,:); 

end 

%Sums  the  fading  paths 
y=c*xd; 

% 

%Adds  in  the  First  path  without  fading 

y=y+xo; 

y=y(l:F*Nt); 

y=reshape(y,Nt,F).'; 
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%CMV2FA 

% 

% 

%Title  :  Complex  Frequency  Array  Generator 
%Author  :  Dave  Roderick,  Naval  Postgraduate  School 
%Modified  By  :  Ahmet  Yasin  ERDOGAN,  Naval  Postgraduate  School 

% 

%Changes  complex  modulation  values  to  frequency  array 

% 

function  X=cmv2fa(N,M) 

[m  n]=size(M); 
if  rem(n,2)==0; 

M=M; 

else 

% 

M=[zeros(m,  1 )  M] ; 
end 

[m  n]=size(M); 

K=round(n/2); 

% 

%Generate  a  matrix  of  zeros  with  m  row  and  N  columns. 
X=zeros(m,N); 

% 

X(:,1:K)=M(:,K+1:2*K); 

X(:,N-K+1:N)=M(:,1:K); 
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%CNV_ENCD 

% 

% 

%Title  :  Convolutional  Encoding 

%Reference  :  Contemporary  Communication  System  using  Matlab 
%John  G.  Proakis  &  Masoud  Salehi. 

% 

function  ce_output=cnv_encd(ce_g,ce_kO,ce_input) 

%  cnv_encd(ce_g,ce_kO,ce_input) 

%  determines  the  output  sequence  of  a  binary  convolutional  encoder 
%  ce_g  is  the  generator  matrix  of  the  convolutional  code 
%  with  ce_nO  rows  and  ce_l*ce_kO  columns.  Its  rows  are  ce_gl,ce_g2,....ce_gn. 
%  ce_kO  is  the  number  of  bits  entering  the  encoder  at  each  clock  cycle. 

%  check  to  see  if  extra  zero  padding  is  necessary 
if  rem(length(ce_input),ce_kO)>0 

ce_input=[ce_input,zeros(size(  1  :ce_kO-rem(length(ce_input),ce_kO)))] ; 
end 

ce_n=length(ce_input)/ce_kO; 

%check  the  size  of  matrix  ce_g 
if  rem(size(ce_g,2),ce_k0)>0 
error('Error,  ce_g  is  not  of  the  right  size.') 
end 

%  determine  ce_l  and  ce_nO 
ce_l=size(ce_g,2)/ce_k0; 

%disp(['The  value  of  ce_l  is:',int2str(ce_l)]); 
ce_nO=size(ce_g,  1 ) ; 

%disp(' ') 

%disp(['The  value  of  ce_nO  is:',int2str(ce_n0)]); 

%add  extra  zeros 

ce_u=[zeros(size(  1  :(ce_l- 1  )*ce_kO)),ce_input,zeros(size(  1  :(ce_l- 1  )*ce_kO))] ; 

%generate  ce_uu,  a  matrix  whose  column  are  the  contents  of 

%conv.  encoder  at  various  cycles. 

ce_ul=ce_u(ce_l*ce_kO:- 1:1); 

for  ce_i=  1  :ce_n+ce_l-  2 

ce_ul=[ce_ul,ce_u((ce_i+ce_l)*ce_kO:-l:ce_i*ce_kO+l)]; 

end 

ce_uu=reshape(ce_u  1  ,ce_l*ce_kO,ce_n+ce_l- 1 ) ; 

%determine  the  ce_output 

ce_output=reshape(rem(ce_g*ce_uu,2),l,ce_n0*(ce_l+ce_n-l)); 
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clc, close  all 
%COFDMSIM 

% 

% 

%Title  :  Simulation  Of  COFDM 

%Author  :  Dave  Roderick,  Naval  Postgraduate  School 

%Modified  By  :  Ahmet  Yasin  ERDOGAN,  Naval  Postgraduate  School 

% 

% 

disp('_ _ 

disp('This  batch  m-file  mns  COFDM  simulations  using  different  channel  models.') 
%fort=input('To  run  the  frequency  version,  enter  l(one),  To  mn  the  time  version,  enter 
O(zero),  or  to  mn  both  enter  2(two):'); 
fort=l;  %frequency  version 

%freqno=input('Enter  the  #  of  OFDM  frequencies  (note  :  must  be  even):'); 
freqno=48; 

%%%%N=input('Enter  the  number  of  EFT  points  (Note  :  This  number  must  be  larger 
than  #  of  OFDM  frequencies):'); 

N=64; 

chnmdl=input('Choose  the  channel  model;  0-(Noise  Free),  l-(AWGN),  21-(Mobile 
Channel- 1),  22-(MC-2),  23-(MC-3),  24-(MC-4)?  :'); 
if  chnmdl=0 

disp('Code  Check  simulation .'); 
sigs=0;loss=0;dop=0;dly=0; 
elseif  chnmdl=l 

disp('AWGN  Channel  simulation .'); 

sigs=input('Enter  the  noise  variance(sigma)  range  or  single  value.  (Ex  lin- 
space(0,0.02,20)or  .003):'); 

loss=0;dop=0;dly=0; 
elseif  chnmdl=21 

disp('Mobile  Channel- 1.'); 

sigs=input('Enter  the  noise  variance(sigma)  range  or  single  value.  (Ex  lin- 
space(0,0.02,20)or  .003):'); 

loss=[0,2. 17,4.34, 6.51,8.69, 10.86, 13.03, 15.20, 17.37,19.54,21.71, 23. 89,26.06,28.23,30.4, 
32.57,34.74,36.92]; 

dop=[15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15]; 

dly=[0,0.05,0. 10,0. 15,0.20,0. 25,0. 30,0. 35,0.40,0.45,0. 50,0. 55,0. 60,0. 65,0.70,0.75,0. 80,0. 
85]; 

elseif  chnmdl=22 

disp('Mobile  Channel-2.'); 

sigs=input('Enter  the  noise  variance(sigma)  range  or  single  value.  (Ex  lin- 
space(0,0.02,20)or  .003):'); 
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loss=[0, 2.17, 4.34, 6.51,8.69, 10.86, 13.03, 15.20, 17.37, 19.54, 21.71, 23. 89,26.06,28.23, 30.4, 
32.57,34.74,36.92]; 

dop=[5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5]; 

dly=[0,0.05,0. 10,0. 15,0. 20,0. 25,0. 30,0. 35,0. 40,0. 45 ,0.50,0. 55,0. 60,0. 65,0. 70,0. 75,0. 80,0. 
85]; 

elseif  chnmdl=23 
disp('Mobile  Channel-3.'); 

sigs=input('Enter  the  noise  variance(sigma)  range  or  single  value.  (Ex  lin- 
space(0,0.02,20)or  .003):'); 

loss=0;dop=0;dly=0; 
elseif  chnmdl=24 

disp('Mobile  Channel-4.'); 

sigs=input('Enter  the  noise  variance! sigma)  range  (Ex:  linspace(0,0.02,20)or  .003):'); 
loss=0;dop=0;dly=0; 
end 

%%%%  Choosing  Interleaver  %%%% 

%allcase=input('Simulate  all  interleaver  cases  (yes)  or  specific  ones(no)?  (l=yes,0=no):'); 

allcase=0; 

if  allcase=l 

disp('All  cases, (0-8), will  be  tested.'); 
cases=[0:8]; 
elseif  allcase==0 

%%%%cases=input('Enter  specific  case  numbers  from  (0  to  8)(Ex  [0  4  5  8]):'); 

cases=0;  %  Block  interleaving  chosen 
end 


if  fort~=2 

if  length(cases)~=l 

casey_n=input('Do  you  want  to  find  optimal  interleaver  case(s)  ?  (l=yes,  0=no):'); 
else 

casey_n=0; 

end 

end 

%%%%% 

totsym=input('Enter  the  total  minimum  number  of  symbols  to  simulate  (Ex  10000):'); 
ro  wno=ceil(totsy  m/freqno) ; 
if  totsym~=(rowno*freqno) 

disp(['Note:The  actual  total  number  of  symbol  to  be  simulated  will  be 
:',int2str(rowno*freqno)]); 
end 

%pry_n=input('For  the  interleaver,  do  you  want  to  calculate  all  possible  intermediate  ma¬ 
trix  dimension  pairs ?(l=yes,0=no):'); 
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pry_n=  1 ; 
pairl=l; 

pair2=rowno*freqno; 
if  pry_n==l 

Intrlvr_pairs=intlvprs(rowno,ffeqno); 
intlvrprs=Intrl  vr_pairs ; 

%%%%  Desired  interleaver  pairs  can  be  entered  %%%% 

%disp(") 

disp('For  these  input  parameters,  all  possible  interleaver  dimension  pairs  are: ') 
disp(Intrlvr_pairs) 
end 

pairs=input([Desired  interleaver  pair?  (Ex  [row  #  col  #]  =  [20  50]  (Note:  enter- 
ing[',int2str(pairl),'  ',int2str(pair2),'],or  [',int2str(pair2),'  ',int2str(pairl),'],  offers  no  inter¬ 
leaving  functionality):']); 

%rintlv=intlvrprs(8,l);cintlv=intlvrprs(8,2) 

rintlv=pairs(l); 

cindv=pairs(2); 

%mary=input('Enter  the  number  of  M- ary  bits,  q  (i.e.  for  256- ary,  q=8):'); 
mary=l; 

%nary=input('Enter  the  number  of  N-ary  bits,q(i.e.  for  16- ary,  q=4):'); 
nary =2;  %  QPSK  is  chosen 
%freqspace=round(  1 6600000/freqno) ; 

freqspace=3 12500;  %  According  to  the  standard  0.3125  Mhz  of  Frequency  spacing. 
%Ng=input('Enter  the  guard  interval  length  (Number  of  sample  points):'); 

Ng=16; 

%ecc=input('Do  you  want  to  include  Reed  Solomon  error  correction  coding  ?  (l=yes, 
0=no):'); 

ecc=0;%  Reed  Solomon  Error  correction  coding  is  not  chosen 
if  ecc==l 

%code=input('Enter  n,k  and  error  correction  block  length  (Ex  [240  200  240]):'); 
code=[240  200  240]; 
n=code(l); 
k=code(2); 
blklgth=code(3); 
elseif  ecc==0 
n=freqno; 
k=freqno; 
blklgth=freqno; 
end 

svals=input('Enter  specific  seed  values,  or  0  for  a  random  seed  (ex  [103  22,  60]  or  [0]):'); 

wait=3;pic=0; 

svect=[]; 

for  mn=l:length(svals); 
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errvect=[]; 
errcase=[]; 
errtot=[]; 
if  min(svals)=0 

rand('seed',sum(100*clock)); 
s=round(abs(rand(  1  )*pi*  1 0*(pic+ 1  )*run)); 
elseif  min(svals)~=0 
s=svals(run); 
end 

svect=[svect,s]; 
for  1=1  iength(cases); 

dispC. _ _') 

disp(['Run#:',int2str(run)]); 
disp(['Seed=',int2str(s)]); 
disp(  ['Interleaver  case=',int2str(cases(l))]); 
if  fort<=l 


chancdl(chnmdl, wait, pic, cases(l),s,freqno,rintlv,cintlv,N,mary, nary  ,n,k,blklgth,Ng,sigs,lo 
ss,dly,dop,freqspace,fort); 
elseif  fort==2 

disp('Frequency  differential  encoding/decoding  simulation...') 

disp('. _ _') 


chancdl(chnmdl, wait, pic, cases(l),s,freqno,rintlv,cintlv,N,mary, nary  ,n,k,blklgth,Ng,sigs,lo 
ss,dly,dop,freqspace,  1); 

disp( ************************************************************) 

disp('Time  differential  encodi ng/dccodi ng  simulation....') 

disp(' _ .’) 


chancdl(chnmdl, wait, pic+12,cases(l),s,freqno,rintlv,cintlv,N,mary, nary  ,n,k,blklgth,Ng,sig 
s,loss,dly,dop,freqspaee,0); 
end 
end 
end 

)  P( 

disp(");disp('Simulation  finished !') 

Seed=svect; 
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%CVDD 

% 

% 

%Title  :  Continuous  Variable  digital  delay  element. 

%Reference  :  C.W.  Farrow, "  A  Continuously  Variable  Digital  Element",  IEEE 
%Intemational  Symposium  on  Circuits  &  Systems, pp.2641-2645, 1988. 

% 

% 

function  [y]=cvdd(x, alpha) 
if  ((nargin~=2)l(nargout~=l)) 
error('ERROR:usage:  y=y=cvdd(x, alpha) ;') ; 
return; 
end 

if  (size(x)~=size(alpha)) 

error('ERROR:x  and  alpha  must  be  the  same  size'); 

return; 

end 

if  (abs(alpha)>0.5) 

errorf'ERROR: alpha  must  be  within  -0.5  and  0.5'); 

return; 

end 

% 

% - - - 

%  Initialization 

% - - - 

% 

%  Initialize  FIR  filter  coefficients  are  in  [1]  (0,0.328  pass  band) 

C0=[-0. 013824  0.054062  -0.157959  0.616394  0.616394  -0.157959  0.054062  -0.013824]; 
Cl=[0.003143  -0.019287  0.1008  -1.226364  1.226364  -0.1008  0.019287  -0.003143]; 
C2=[0.055298  -0.216248  0.631836  -0.465576  -0.465576  0.631836  -0.216248  0.055298]; 
C3=[-0. 012573  0.077148  -0.403198  0.905457  -0.905457  0.403198  -0.077148  0.012573]; 

% 

% 

%  4  parallel  FIR  and  add  together  based  on  [1] 

% . . . . . . 

y0=filter(C0,  [  1  ]  ,x) ; 
yl=filter(Cl,[l],x); 
y2=filter(C2,[l],x); 
y3=filter(C3,  [  1  ]  ,x) ; 

% 

y=alpha.*y3; 

y=alpha.*(y+y2); 

y=alpha.*(y+yl); 

y=y+y0; 
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%DECDRCDL 

% 

% 

%Title  :  COFDM  Decoder  With  Deinterleaveing 
%Author  :  Dave  Roderick,  Naval  Postgraduate  School 
%Modified  By  :  Ahmet  Yasin  ERDOGAN,  Naval  Postgraduate  School 

% 

%-— 

func- 

tion[outmsg,viterbi_output_bit,random_msg,random_bit,M,MM]=decdrcdl(pic,dcase,K,F 

a,nsymno.freqno,rdintlv,cdintlv.mary.nai'y,foit.B_random) 

%  ’ 

M=fa2cma(K,Fa) ; 

Cmplx_mod_vals=M; 

% 

naryp=nary; 

[s,MM]=dfdcdrft(naryp, nary  ,M, fort); 

[L,cc]=size(s); 
strans=s'; 
svect=strans(:).'; 
corrs=s  vect(  1 :  nsy  mno) ; 

% 

nsymno; 

Br=bm(mary,mb(nary  ,corrs)) ; 
lengthBi -length(Br) ; 
mi  ndr=rcm(  length  ( Br),ffeq  no) ; 
if  rmndi^=0; 

Bi-Br; 

elseif  rmndr~=0; 

Bi-Bif  1  :(lengthBr-  rmndr)) ; 
end 

rcvd=(reshape(Br,freqno,length(Br)/freqno))'; 

Rcvd_Indv_Ary=rcvd; 

% 

[Br  Bc]=size(rcvd); 

SYNC=[]; 

sr=rcvd'; 

si=sr(:)'; 

sd=cdldlv(rdintlv,cdintlv,dcase,si,SYNC); 
received=reshape(sd,Bc,Br)'; 
viter_G=[l  0  1  1  0  1  1  ;1  1  1  1  0  0  1]; 
viter_k=l; 

binary_value=mb(mary  ,sd) ; 

[viterbi_output,survivor_sta,cumul_metrix]=viterbi(viter_G,viter_k,binary_value); 
mary_dec=bm(mary  ,viterbi_output) ; 
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viterbi_output_bit=viterbi_output; 

%outmsg=reshape(sd,Bc,Br)'; 

% 

random_bit=B_random; 
random_msg=bm(mary  ,random_bit) ; 
[Brow  Bcol]=size(random_msg); 

% 

outmsg=reshape(mary_dec,Bcol,Brow)'; 

Sink_Msg=outmsg; 
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%DECI2BIN 

% 

% 

%Title  :  Decimal  to  Binary  Converter 

%Author  :  Tan  Kok  Chye,  Naval  Postgraduate  School 

% 

function  y=deci2bin(x,l) 

y=zeros(l,l); 

vi=l; 

while  x>=0  &  vi<=l 

y(vi)=rem(x,2); 

x=(x-y(vi))/2; 

vi=vi+l; 

end 

y=y(i:-  i:i); 
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%DFDCDRFT 

% 

% 

%Title  :  Complex  Number  Demodulator  &  Frequency/Time  Differential  Decoder 
%Author  :  Prof.  Paul  H.  Moose,  Naval  Postgraduate  School 
%Modified  by  :  Ahmet  Yasin  ERDOGAN,  Naval  Postgraduate  School 

% 

%-- 

function  [s,M]=dfdcdrft(qp,q,MD,fort) 
if  fort==0  %Time  Differential  decoding 

% 

MD=MD'; 

[m  n]=size(MD); 

% 

%  Perform  a  looping  routine  to  find  the  phase  differences  between  adjacent  values  in  the 
%  array ,MD,and  put  these  calculated  values  into  array ,M. 
for  1=1  :m 
for  j=l:n-l 

M(l,j)=MD(l,j+l)*conj(MD(l,j)); 

end 

end 

% 

%Transpose  the  array  back  to  its  original  form 
M=M'; 

% 

%  Calculate  the  number  of  M-ary  symbols  based  upon  the  exponent  qp,then  use  this 
number 

%  to  find  the  number  of  equally  spaced  phases  in  a  unit  circle. 

N=2Aqp; 

dph=2*pi/N; 

% 

%  Divide  the  phase  arguments  of  elements  in  M,  by  the  equal  phases  generated  by  dph. 
phn=angle(M)  ./dph; 

% 

%  Calculate  the  phase  sector  number  by  finding  the  remainders. 

s=rem(round(phn)+N,N); 

elseif  fort==l  %  Frequency  Differential  decoding 

% 

%  Transpose  the  modulation  array,  and  find  the  dimensions 
[m,n]=size(MD); 

MD=MD(:,2:n); 

[m  n]=size(MD); 

% 

%  Perform  a  looping  routine  to  find  the  phase  differences  between 
%  adjacent  values  in  the  array,  MD,  and  put  these  calculated  values  into  array ,M. 
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for  1=1  :m 
for  j=l:n-l 

M(l,j  )=MD(l,j+ 1  )*conj  (MD(l,j )); 
end 
end 
N=2Aqp; 
dph=2*pi/N; 

% 

%  Calculate  the  phase  sector  number  by  finding  the  remainders. 

phn=angle(M)  ./dph; 

s=rem(round(phn)+N,N); 

end 
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%DIFCDRFT 

% 

% 

%Title  :  Complex  Number  Modulator  &  Frequency/Time  Differential  Encoder 
%Author  :  Prof.  Paul  H.  Moose,  Naval  Postgraduate  School 
%Modified  by  :  Ahmet  Yasin  ERDOGAN,  Naval  Postgraduate  School 

% 

function  MD=difcdrft(q,m,fort) 
if  fort==0  %Time  differential  encoding 
%  M-ary  alphabet  size 
N=2Aq; 

%  Determine  the  number  of  equal  phases  based  upon  the  m-ary  symbol  length 
dph=2*pi/N; 

%  Find  the  size  of  the  input  symbol  matrix 
[rr  n]=size(m); 

% 

%  Perform  the  time  differential  encoding  of  phase  values  by  cumulative  summation, 

%  down  one  column  at  a  time  across  the  entire  matrix.  This  function  generates  a  matrix, 
for  k=l:n 

md=cumsum(  m(  :,k)); 

%  Generate  the  complex  numbers  with  corresponding  phase  values. 

MD(:  ,k)=exp(i*dph.*md) ; 
end 

% 

%  Inject  the  reference  row  of  ones  (zero  phase)  at  top  of  output  matrix  for 
%  differential  encoding  synchronization 
MD=[ones(l,n);  MD]; 

elseif  fort==l  %  Frequency  Differential  encoding 

% 

%  M-ary  alphabet  size 

N=2Aq; 

dph=2*pi/N; 

%  Find  the  size  of  the  input  symbol  matrix 
[rr  n]=size(m); 

% 

md=cumsum(m'); 

md=md'; 

% 

%  Generate  the  complex  numbers  with  corresponding  phase  values. 

MD=exp(i*dph.  *md) ; 

% 

%  Inject  the  reference  row  of  ones  (zero  phase)  at  top  of  output  matrix  for 
%  differential  encoding  synchronization. 

MD=[ones(rr,2)  MD]; 
end 
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%DIFFCHKR 

% 

% 

%Title  :  Differential  Encoder/Decoder  Checker 
%Author  :  Dave  Roderick,  Naval  Postgraduate  School 
%Modified  By  :  Ahmet  Yasin  ERDOGAN,  Naval  Postgraduate  School 

% 

%— - 

function  diffchkr(s,symno,freqno,mary,nary) 

fort=input(Eor  the  frequency  version,  enter  1  (one);  for  the  time  version,  enter  0  (zero):'); 

% 

B=marymsg(mary,s,symno,freqno); 

Rndm_m_ary_msg=B ; 

% 

ml=bm(nary,mb(mary,B)); 

lengthml=length(ml); 

m=(reshape(ml,lengthml/symno,symno))'; 

N_ary_msg=m; 

% 

if  fort==l 

disp(");disp('Frequency  Differential  Encoding/Decoding  version') 

% 

%Freq.  Diff.  Enc. 

% 

MDD=difcdrf(mary,m); 
elseif  fort~=l 

disp(");disp('Time  Differential  Encoding/Decoding  version') 

% 

MDD=difcdrt(mary  ,m) ; 
end 

% 

maryq=mary; 
if  fort==l 

% 

[s  M]=difdcdrf(maryq,mary,MDD); 
elseif  fort  ~=1 

% 

[s  M]=difdcdrt(maryq,mary,MDD); 
end 

% 

%Check  results  for  correctness. 

[error_no,freqerrs,errmx,rowerrs]=check(0,m,s,freqno,freqno,freqno); 
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%DLINE 

% 

% 

%Title  :  UHF  Channel  Delay  Line  Generator 

%Author  :  Prof.  Paul  H.  Moose,  Naval  Postgraduate  School 

% 

function  xd=dline(x,d) 
x=x.'; 

dmax=max(d); 
dmin=min(d); 
nmin=floor(drnin); 
nmax=ceil(dmax) ; 
x=[x;zeros(nmax+3, 1)] ; 

N=length(x); 

Nd=length(d); 

% 

for  n=l:Nd; 
di=d(n); 

D=floor(di); 

dcld=di-D; 

xd(:,n)=cvdd(x,deld-  .5); 
xd(:,n)=[zeros(D,l);xd(l  :N -D,n)]; 
end 

xd=xd.'; 

[rr,cc]=size(xd); 
xd=xd(:  ,4+nmin:cc) ; 
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%FA2CMA 

% 

% 

%Title  :  Frequency  Array  To  Complex  Modulation  Array  Converter 
%Author  :  Prof.  Paul  H.  Moose,  Naval  Postgraduate  School 
%Modified  by  :  Ahmet  Yasin  ERDOGAN,  Naval  Postgraduate  School 

% 

%— - 

function  Mm=fa2cma(K,X) 

[m  n]=size(X); 

Mm(:,l:K)=X(:,n-K+l:n); 

Mm(:,K+l:2*K)=X(:,l:K); 

Cmplx_mod_vals=Mm; 
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%FROFST 

% 

%Title  :  Frequency  Offset 

%Author  :  Ahmet  Yasin  ERDOGAN,  Naval  Postgraduate  School 

% 

% 

function  [xofstd]=frofst(x,y) 

[sati  sutu]=size(x); 
x=x.'; 
x=x(:).'; 
wo=length(x); 
wu=l:wo; 
yd=wu*y 

xofst=x.  *exp(i*(2 .5  *pi) .  *yd) ; 
xofst=xofst(  1 :  sutu*  sati) ; 
xofst=reshape(xofst,sutu,sati).'; 
xofstd=xofst; 
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%INTLVCHK 

% 

% 

%Title  :  Interleaver/Deinterleaver  Verifier 

%Author  :  Dave  Roderick,  Naval  Postgraduate  School 

%Modified  By  :  Ahmet  Yasin  ERDOGAN,  Naval  Postgraduate  School 

% 

%— - 

function  intlvchk(s,symno,lfeqnojintlv,cintlv,mary,dcase) 
multiples=mltpl(symno,freqno); 

Intrlvr_nbr_mltpls=multiples; 

% 

if  (symno*freqno)~=(rintlv*cintlv) 

disp('ERROR:  The  interleaver  matrix  dimensions  are  not  compatible  with  the  message 
array  size.') 

disp('  Possible  matrix  dimensions  are:');disp(") 
dispfmultiples) 

disp('Note:  The  selected  matrix  dimensions  can  not  accomodate  the  message  array.') 
disp('  In  this  case  the  number  of  rows  times  the  number  of  columns  is:');disp(") 
disp(symno*freqno) 
elseif(symno*fieqno)/(rintlv*cintlv)=  1 

% 

B=maiymsg(mary,s,symno,freqno) ; 

Random_msg=B 

% 

SYNC=[]; 

[Br  Bc]=size(B); 

Bt=B'; 

Bvect=Bt(:)'; 

si=cdhlv(rintlv,cintlv,dcase,Bvect,SYNC); 

Bi=reshape(si,Bc,Br)'; 

Interleaved_array=Bi 

% 

[Br  Bc]=size(Bi); 

SYNC=[]; 

sr=Bi'; 

si=sr(:)'; 

sd=cdldlv(rintlv,cintlv,dcase,si,SYNC); 

Bd=reshape(sd,Bc,Br)'; 

Deinterleaved_array=Bd 

% 

[error_no,fi'eqerrs,errmx,rowerrs]=check(0,B,Bd,freqno,freqno,freqno); 

end 
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%INTLVPRS 

% 

% 

%Title  :  Intermediate  Matrix  Interleaver  Dimension  Pairs 
%Author  :  Dave  Roderick,  Naval  Postgraduate  School 
%Modified  By  :  Ahmet  Yasin  ERDOGAN,  Naval  Postgraduate  School 

% 

%— - - - - - - - - - 

function  pairs=intlvprs(n,m) 
prod=n*m; 
multvect=[l]; 
for  i=2:prod; 

remdr=rem(prod,i) ; 
if  remdi-=0 

multvect=[multvect  ij; 
else 

multvect=multvect; 

end 

mult=multvect; 

end 

lngth=length(mult); 
nhi-mulUlngth); 
result=[l  nbr]; 
for  i=2:lngth; 

cmtpr=[mult(i)  nbr/mult(i)]; 
result=[result;cmtpr] ; 
end 

pairs=result; 
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%ITDA 

% 

% 

%Title  :  Frequency  Domain  Samples  Without  Guard  Interval 
%Author  :  Prof.  Paul  H.  Moose,  Naval  Postgraduate  School 

% 

% 

function  Y=itda(Ng,y) 

% 

[L  Nt]=size(y); 

%  Remove  the  guard  interval  for  channel  compensation. 

% 

y=y(:,Ng+l:Nt); 

%  Take  the  FFT  of  array,  y 

% 

Y=fft(y.').'; 
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%JAKES 

% 

% 

%Title  :  Multi-path 

%Author  :  Prof.  Roberto  Cristi,  Naval  Postgraduate  School 
%Modified  by  :  Ahmet  Yasin  ERDOGAN,  Naval  Postgraduate  School 

% 

%-— - - - - - 

function  g=jakes(Fd,Fs,N); 

%  g=jakes(Fd,Fs,N) 

%  generate  a  random  vector  of  length  N  with  jakes  spectrum 
%  Fd=doppler  frequency 
%  Fs=sampling  frequency  (Fs>10*Fd  required) 

%  N=vector  length 

%  g=generated  random  vector 

%  The  vector  g  has  unit  power,  ie.  g'*g*Fs/N=l 


F=round(log2(Fs/(  10*Fd))); 

F0=Fs/(2AF);  %  samphng  frequency  for  generating  the  random  sequence 

% 

N0=2*ceil(N/(2AF)); 

Nf=5 12;  %  number  of  frequency  components  in  spectrum  definition 

Nh=256;  %  length  of  FER  filter 

%  FIR  Filter  impulse  response 

w=0.54-0.46*cos(2*pi*(0:Nh-  1)/Nh);  %  hamming  window 

fd=Fd/F0;  %  digital  doppler  frequency 

f=0:(l/Nf):l-(l/Nf);  %  vector  of  digital  frequencies  (l/2=Nyquist  Freq.) 

kd=floor(fd*Nf);  %  index  for  doppler  frequency  fd 

H0(  1  :kd)=sqrt(  1  -  (f(  1  :kd)/fd) .  A2) ; 

I=find(H0==0);  H0(I)=0.0001*ones(size(I)); 

H=zeros(l,Nf); 

H(  1  :length(H0))=(  1  ./HO) ; 

H=H.*exp(-j*2*pi*f*Nh/2); 
hO=real(ifft(H)); 
h=hO(  1  :length(w)).*w; 
end 

%  Generate  time  varying  taps  at  low  sampling  freq  Fs_ch 
seed=14; 

randn('state',seed); 
a=randn(  1  ,N0+Nh) ; 
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seed=seed+6; 
b=randn(  1  ,NO+Nh) ; 

%x=randn(  1  ,NO+Nh)+j  *randn(  1  ,NO+Nh) ; 
x=a+j*b; 

g=filter(h,l,x); 

g=g(Nh+l:Nh+NO);  %  steady  state  response 


%  Upsample  to  Fs=(2AL)*F0  in  L  stages 

for  m=l:L 

F0=2*F0; 

g=reshape([g;zeros(size(g))],  l,2*length(g)); 
omegad=2*pi*Fd/F0; 

Domega=(pi/2)-  omegad; 
omegac=((pi/2)+omegad)/2 ; 

M=ceil(((8  *pi/Domega)- 1  )/2) ; 

Nfilt=2*M+ 1 ;  %  filter  order  using  hamming  window 

nt=0:Nfilt- 1 ; 

w=0.54-0.46*cos(2*pi*nt/Nfilt);  %  hamming  window 

hm(  1  :M)=sin(omegac*(nt(  1  :M)-  M))  ./(pi*(nt(  1  :M)-  M)) ; 
hm(M+ 1  )=omegac/pi ; 

hm(M+2:Nfilt)=sin(omegac*(nt(M+2:Nfilt)-M))./(pi*(nt(M+2:Nfilt)-M)); 

hm=hm(l  :Nfilt).*w; 

g=filter(hm,l,g); 

g=g(Nfilt:length(g));  %  steady  state 
end 

%  Normalize 
g=g(l:N); 

K=g*g'/N; 

g=g/sqrt(K); 
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%MARYMSG 

% 

% 

%Title  :  M-ary  Message  Test  Pattern  Generator 
%Author  :  Tan  Kok  Chye,  Naval  Postgraduate  School 
%Modified  by  :  Ahmet  Yasin  ERDOGAN,  Naval  Postgraduate  School 

% 

%— - 

function  [vmary_ce4'andom_bit]=marymsg(q,n,m) 

% 

[random_bit]=msg(n*m*q*.5-6); 

conv_g=[l  0  1  101  1;1  1  1  1001]; 
conv_k0=l; 

conv_output=cnv_encd(conv_g,conv_kO,random_bit); 
decml=bm(q,conv_output) ; 
deciml=decml(  1,1:  (n*m)) ; 
vmary_ce=(reshape(deciml,m,n))' ; 
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%MB 

%- . . . - . — - 

% 

%Title  :  M-ary  To  Binary  Converter 

%Author  :  Prof.  Paul  H.  Moose,  Naval  Postgraduate  School 

% 

function  [b]=mb(q,m) 

% 

row=size(m,l); 

col=size(m,2); 

% 

m=reshape(m',  1  ,(row*col)) ; 

% 

b0=rem(m,2); 

m=(m-b0)./2; 

B=bO; 

% 

for  j=l:q-l 

bj=rem(m,2); 

m=(m-bj)./2; 

B=[B;bj]; 

end 

% 

b=B(:)'; 

binary=b; 
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%METRIC 

% 

% 

%Title  :  Viterbi  Hard  Decision  Decoding  metric 
%Author  :  Tan  Kok  Chye,  Naval  Postgraduate  School 

% 

function  distance=metric(v_x,v_y) 

if  v_x==v_y 

distance=0; 

else 

distance=l; 

end 
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%MLTPL 

%- . — . . . 

% 

%Title  :  Common  Multiples 

%Author  :  Tan  Kok  Chye,  Naval  Postgraduate  School 

% 

function  [mult]=mltpl(n,m) 

max=n*m; 

multvect=[l]; 

% 

for  i=2:max; 
rcmdi-rcm(maxj) ; 
if  remdi-=() 
multvect=[multvect  ij; 
else 

multvect=multvect; 

end 

mult=multvect; 

end 
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%MSG 

% 

% 

%Title  :  Message  Test  Pattern  Generator 

%Author  :  Prof.  Paul  H.  Moose,  Naval  Postgraduate  School 

%Modified  by  :  Ahmet  Yasin  ERDOGAN,  Naval  Postgraduate  School 

% 

%— - 

function  u=msg(k) 

% 

%rand('uniform') ; 

% 

%temp=rand('seed') ; 

% 

%rand('seed',s); 

% 

u=randint(l,k); 
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%NXT_STAT 

% 

% 

%Title  :  Next  State 

%Author  :  Tan  Kok  Chye,  Naval  Postgraduate  School 

% 

function  [next_state,memory_contents]=nxt_stat(current_state, input, v_L,v_k) 
binary_state=deci2bin(current_state,v_k*(v_L- 1 )) ; 
binary_input=deci2bin(input,v_k) ; 

next_state_binary=[binary_input,binary_state(l:(v_L-2)*v_k)]; 

next_state=bin2deci(next_state_binary); 

memory_contents=[binary_input,binary_state]; 
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%OFST 

% 

% 

%Title  :  Channel  Offset 

%Author  :  Prof.  Paul  H.  Moose,  Naval  Postgraduate  School 
%Modified  by  :  Ahmet  Yasin  ERDOGAN,  Naval  Postgraduate  School 

% 

%— - 

function  xo=ofst(e,N,x) 

[m  Nt]=size(x); 

xo=x.'; 

x=x.'; 

x=x(:); 

x=x.'; 

Nt=length(x); 

1=1  :Nt; 

%Creating  the  offset  frequency 

% 

ex=x.*exp(i*(2*pi/N)*e.*l); 

xo(:)=x; 

xo=xo.'; 
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%phsns2 

% 

% 

%Title  :  Phase  Noise 

%Author  :  Ahmet  Yasin  ERDOGAN,  Naval  Postgraduate  School 

% 

%This  m.file  creates  phase  noise  as  a  Wiener  Process  for  all  time  samples 

% 

function  phsnsig=phsns2(xmtifft,freqspace,perctB) 
perctB; 

[str  stn]=size(xmtiffl); 

T  =str*stn; 

B=freqspace*perctB ; 
randn(' state',  100)  ; 
dW  =  zeros(l,stn*str); 

W  =  zeros(l,stn*str); 

% 

dW(l)  =  normrnd(0,sqrt(2*B)); 

W(l)  =  dW(l); 
for  j  =  2:(stn*str) 
dW(j)  =  normrnd(0,sqrt(2*B*j)); 

W(j)  =  (W(j-1)  +  dW(j)); 
end 

xmtifft=xmtifft.'; 

xmtifft=xmtifft(:).'; 

xmtifftpn=xmtifft.*exp(i*(2*pi*4*  1 0A(-  6)).  *W) ; 

xmtifftpn=reshape(xmtifftpn,stn,str).'; 

phsnsig=xmtifftpn; 

%  phase  noise  plot 
figure(  14) 

plot([0:T],[0,W],'r-') 

xlabel('t') 

ylabel('PN(t)') 

title('Phase  Noise  for  all  Time  Domain  Samples') 
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%phsns 

% 

%Title  :  Phase  Noise 

%Author  :  Ahmet  Yasin  ERDOGAN,  Naval  Postgraduate  School 

% 

%This  m.file  creates  phase  noise  as  a  Wiener  Process  for  each  OFDM  symbol 

% 

function  phsnsig=phsns(xmtifft,freqspace,perctB) 

[str  stn]=size(xmtifft); 

T  =stn; 

B=freqspace*perctB ; 

for  i=l:str 

randn('state',100)  ; 
dW  =  zeros(l,stn); 

W  =  zeros(l,stn); 

dW(l)  =normmd(0,sqrt(2*B)); 

W(l)  =  dW(l); 
for  j  =  2:stn 

dW(j)  =  normrnd(0,sqrt(2*B*j)); 

W(j)  =  (W(j- 1)  +  dW(j)); 
end 

xmtifftpn(i,:)=xmtifft(i,:).*exp(i*(2*pi*4*10A(-6)).*W); 

end 

phsnsig=xmtilftpn; 

%  phase  noise  plot 
figure(  14) 

plot([0:T],[0,W],'r-') 

xlabel('t') 

ylabel('PN(t)') 

title('Phase  Noise  Created  for  Each  OFDM  Symbol') 
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%RAY_DOP 

% 

% 

%Title  :  Rayleigh  Doppler 

%Author  :  Prof.  Paul  H.  Moose,  Naval  Postgraduate  School 
%Modified  by  :  Ahmet  Yasin  ERDOGAN,  Naval  Postgraduate  School 

% 

%— - - - - — - . . . 

function  c=ray_dop(s,M,N,es) 

m=0:M-l; 

randn('seed',s+10); 

prl=randn(l,20); 

randn('seed',s+20); 

pim=i*randn(  1 ,20) ; 

p=prl+pim; 

p=p/(40A.5); 

rand('seed',s+30); 

e=rand(l,20); 

e=es*cos(2*pi*(e-  .5)); 

E=exp(i*2*pi*e'*m/N); 

c=p*E; 
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%ROTM 

% 

% 

%Title  :  Rotate  Vector 

%Author  :  Prof.  Paul  H.  Moose,  Naval  Postgraduate  School 

% 

function  [vp,vn]=rotm(v,m) 

L=length(v); 

m=rem(m,L); 

ii=(l:L)-l; 

isp=rem(ii-  m+L,L)+ 1 ; 
isn=rem(ii+m+L,L)+l ; 
vp=v(isp); 
vn=v(isn); 
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%TDA 

% 

% 

%Title  :  lime  Domain  Samples  With  Guard  Interval  Precursor 
%Author  :  Prof.  Paul  H.  Moose,  Naval  Postgraduate  School 
%Modified  by  :  Ahmet  Yasin  ERDOGAN,  Naval  Postgraduate  School 

% 

%— - 

function  x=tda(Ng,X) 

[m  N]=size(X); 

% 

%  Perform  inverse  EFT  on  frequency  values  in  array ,X 
x=ifft(X.'); 

%  Add  precursor  of  Ng  samples  to  the  beginning  of  the  time  domain  array  for  channel 
%  compensation. 

% 

x=x.'; 
if  Ng=0 
x=x; 
else 

x=[x(:,N-Ng+l:N)  x]; 
end 


122 


%VITERBI 

% 

% 

%Title  :  Viterbi  Decoder 

%Reference  :  Contemporary  Communication  System  using  MatLab 
%John  G.  Proakis  &  Masoud  Salehi. 

% 

func- 

tion[decoder_output,survivor_state,cumulated_metric]=viterbi(v_G,v_k,channel_output) 

% 

v_n=size(v_G,l); 

%  check  the  sizes 
if  rem(size(v_G,2),v_k)~=0 
error('Size  of  v_G  and  v_k  do  not  agree') 
end 
% 

if  rem(size(channel_output,2),v_n)~=0 
error('channel  output  not  of  the  right  size') 
end 

v_L=size(v_G,2)/v_k; 
number_of_states=2  A((v_L- 1 )  *  v_k) ; 

%generate  state  transition  matrix,  output  matrix,  and  input  matrix 
for  vj=0:number_of_states- 1 
for  v_l=0:2Av_k- 1 

[next_state,memory_contents]  =nxt_stat(v  _j  ,v_l,v_L,v_k) ; 
input(v  j + 1  ,next_state+ 1  )=v_l; 
branch_output=rem(memory_contents*v_G',2); 
nextstate(v  J+ 1 ,  v_l+ 1  )=next_state ; 
output(v  j + 1  ,v_l+ 1  )=bin2deci(branch_output) ; 
end 
end 

state_metric=zeros(number_of_states,2); 

depth_of_trelhs=length(channel_output)/v_n; 

channel_output_matrix=reshape(charmel_output,v_n,depth_of_trellis); 

survivor_state=zeros(number_of_states,depth_of_trelhs+l); 

%start  decoding  of  non- tail  channel  outputs 
for  v_i=  1  :depth_of_trellis  -  v_L+ 1 
flag=zeros(  1  ,number_of_states) ; 
if  v_i<=v_L 

step=2A((v_L-  v_i)*v_k); 

else 

step=l; 

end 

for  vj=0:step:number_of_states- 1 
for  v_l=0:2Av_k- 1 
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branch_metric=0; 

binary_output=deci2bin(output(v  _j + 1  ,v_l+ 1 ),  v_n) ; 
for  v_ll=  1  :v_n 

branch_metric=branch_metric+metric(channel_output_matrix(v_ll,v_i),binary_output(v_ 

m 

end 

if((state_metric(nextstate(v  _j+ 1  ,v_l+ 1 )+ 1 ,2)>state_metric(v  J+ 1 , 1  )+branch_metric)lflag( 
nextstate(v  J + 1 ,  v_l+ 1 )+ 1  )==0) 

state_metric(nextstate(v  _j + 1  ,v_l+ 1 )+ 1 ,2)=state_metric(v  J +1,1  )+branch_metric ; 
survivor_state(nextstate(  v  J + 1  ,v_l+ 1 )+ 1  ,v_i+ 1  )=v J ; 
flag(nextstate(v  _j+ 1  ,v_l+ 1 )+ 1  )=  1 ; 
end 
end 
end 

state_metric=state_metric(:,2:- 1:1); 
end 

%start  decoding  of  the  tail  channel_outputs 
for  v_i=depth_of_trellis  -  v_L+2  :depth_of_trellis 
flag=zeros(  1  ,number_of_states) ; 

last_stop=number_of_states/(2A((v_i-depth_of_trellis+v_I^2)*v_k)); 
for  vj=0:last_stop- 1 
braneh_metric=0; 

binary_output=deci2bin(output(v  j + 1 , 1 ),  v_n) ; 
for  v_ll=l:v_n 

branch_metric=branch_metric+metric(channel_output_matrix(v_ll,v_i),binary_output(v_ 

ii)); 

end 

if ((state_metric(nextstate(v  J + 1 , 1 )+ 1 ,2)>state_metric(v  j +1,1  )+branch_metric)  lflag(next 
state(vj+l,l)+l)==0) 

state_metric(nextstate(v  _j+ 1 , 1 )+ 1 ,2)=state_metric(v  _j+ 1 , 1  )+branch_metric; 
survivor_state(nextstate(v  J + 1 , 1 )+ 1  ,v_i+ 1  )=v  J ; 
flag(nextstate(v  J + 1 , 1 )+ 1  )=  1 ; 
end 
end 

state_metric=state_metric(:,2:- 1:1); 
end 

%generate  the  decoder  output  from  the  optimal  path 
state_sequence=zeros(  1  ,depth_of_trellis+ 1 ) ; 

state_sequence(  1  ,depth_of_trellis)=survivor_state(  1  ,depth_of_trelhs+ 1 ) ; 
for  v_i=l:depth_of_trellis 
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state_sequence(  1  ,depth_of Jrellis  - 

v_i+ l)=survivor_state((state_sequence(  1  ,depth_of_trellis+2-  v_i)+l  ),depth_of_trellis  - 

v_i+2); 

end 

decoder_output_matrix=zeros(v_k,depth_of_trellis  -  v_L+ 1 ) ; 
for  v_i=l:depth_of_trellis-v_L+l 

dec_output_deci=input(state_sequence(  1  ,v_i)+ 1 ,  state_sequence(  1 ,  v_i+ 1 )+ 1 ) ; 
dec_output_bin=deci2bin(dec_output_deci,v_k); 
decoder_output_matrix(:,v_i)=dec_output_bin(v_k:- 1: 1)'; 
end 

decoder_output=reshape(decoder_output_matrix,  1  ,v_k*(depth_of_trellis  -  v_L+ 1 )) ; 
cumulated_metric=state_metric(  1,1); 
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